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SUMMARY: The result derived in a previous paper") for the primary effect of 

the jet mixing on the thrust is confirmed, both by an investigation using 

momentum principles and by a quantitative examination of the sink effect. 

Boundary layer behaviour still cannot be predicted with accuracy, but it is 

suggested that, even under cruising conditions, a jet flapped aerofoil may 
exhibit unconventional properties. 


Introduction 


This paper contains a fuller theoretical investigation of certain aspects of the 
paper “ Mixing and the Jet Flap” which appeared in the previous issue of The 
Aeronautical Quarterly. Within the scope of this earlier paper it was not possible 

mto analyse the flow on the basis of momentum, even though this would have been 
the more conventional approach, and neither was it possible, when suggesting a 
picture for the flow, to examine the quantitative arguments for an equivalent sink 
distribution. In the present paper an attempt is made to remedy these defects and 
also to develop the discussion of the mixing with the boundary layer. 


NOTATION 
T. net thrust on a conventional aerofoil 

net thrust on a jet flapped aerofoil 
static pressure in the undisturbed main stream 
static pressure at the position of the mixing 
velocity of the undisturbed main stream 
main stream velocity corresponding to p, 

= U,, jet velocity, at pressure p, before mixing 


jet velocity at pressure p, before mixing 
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ua Velocity difference between the jet and the local main stream 
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density, or densities, of main stream fluid 
density of jet fluid 

mass flow rate of the jet 

sink strength per unit span 


velocity component of the external flow perpendicular to the 
boundary of the jet mixing region 


ratio of the mass flow rates of the main stream and of the jet in the 
mixing process at pressure p, 


distance measured perpendicular to the jet or boundary layer 


local velocity in the jet or boundary layer 


f . (1 - 7 dy, momentum thickness of the boundary layer 
1 1 
0 


co 


(1 ~ - ) dy, displacement thickness of the boundary layer or jet 
1 


chord of the aerofoil 


skin friction 
where p is the local static pressure 


velocity difference between the jet and the main stream at pressure p, 


= Ua,o, Velocity difference between the jet and the main stream at 


pressure p, 


(U,/U,)- 1 


= 1 


“contraction ” of the external flow due to the jet mixing 


refers to conditions at the position of mixing 
refers to conditions in the undisturbed main stream 
refers to the jet. 


refers to values after mixing 


refers to conditions at infinity, which is at pressure p, 
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MIXING AND THE JET FLAP 


2. Momentum Principles Applied to the Mixing Between the Jet 
and the Main Stream 


In this section momentum principles will be used, instead of an energy balance, 
to confirm the result of Ref. 1, equation (15), that 


1) (1 = PU, 
T; -T, U, 


(1+ 


It is still assumed that the flow directions are parallel in the mixing region—a 
fundamental simplification which has already been discussed“’—and in this and 
certain other respects the derivation will be similar to the conventional analysis for 
the mixing of two streams in pipe flow”. 


(1) 


The calculation is based on the following model. The jet emerges into a region 
where the static pressure is p, and there mixes with n times its own mass flow of the 
main stream; the pressure remains constant during the mixing and momentum is 
conserved in the direction of the jet. When this portion of the mixing is complete 
the flow passes without further mixing to the undisturbed pressure p, at infinity, the 
Bernoulli equation of constant energy applying along any streamline. The mixing 
which occurs at infinity is ignored, as it does not affect the thrust on the aerofoil. 
Moreover the whole flow is treated as incompressible, one aspect of this simplifica- 
tion being that temperature changes due to kinetic energy dissipation are neglected; 
it is further assumed that the jet has a uniform or square velocity profile before 
mixing, while the jet plus entrained air has a uniform velocity profile after mixing. 
This again is a simplification which is discussed in Ref. 1. The model just described 
is compared with a conventional aircraft, where effectively all the mixing occurs 
sufficiently far from the wing for it to be ignored when calculating the thrust. For 
each flow it is supposed that the jet air is taken from an intake on the aerofoil. 


The derivation assumes the following theorem : — 


The net thrust on an aerofoil is equal to the net flux of momentum of the 
flow at infinity. 


The proof of the theorem is beyond the scope of this paper and it will be assumed 
that it is acceptable either on intuitive grounds—for the conventional flow its truth 
is apparent—or from a simple “ momentum box ” argument, or as an extension of 
the standard theorem equating aerofoil drag to the momentum deficiency in the 
wake at infinity. Application of the theorem reduces the problem to that of finding 
the net momentum of the jet flap stream at infinity after mixing has occurred at 
pressure p,. 


As in the previous paper, the symbol wu, is used to denote the velocity difference 
between the jet and the main flow, m represents the mass flow of the jet before 


” 


mixing, and (n+1)m the mass flow after mixing. Also suffix “n” is used to 
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indicate that mixing has occurred, while suffix “(n=0)”, referring to flow without 
mixing, is adopted when necessary for the avoidance of ambiguity. In this notation, 
and by the theorem quoted, the thrust T. for the conventional aircraft is 


ic= (n=0) 
while, by the same theorem, the thrust 7; for the jet flap is 
T,=(n + 1) muta, 


The proportional increase of thrust is therefore 


Ua, (n=0) 
The jet velocity before mixing, i.e. Uj, ,. is given by the Bernoulli relationship 
where p,+4pU ,?=p,+4p U,’, 
i.e. U,, 2=U,, 2+ —U2-U,’). 


J 


The jet velocity after mixing is determined by constancy of total momentum in the 
direction of the jet: — 


mU;, ,+nmU,=(n+1) mUj, 


which, with U;, , expressed as (U, + ua, ,), reduces to 


Uas 


Uds 1s n= (n +1) . (4) 
The density of the jet after mixing can be shown to be given by 
gs 


p; 


This result follows from a simple weighted mean when the jet and main stream 
consist of different gases at the same temperature. It similarly follows for streams 
of chemically the same gas at different temperatures, provided that the gas laws 
hold so that volume flow is conserved. 


Knowing the velocity and density of the jet after mixing, the velocity of the 
jet at infinity is given from the Bernoulli relations by the equation corresponding 
to equation (3), i.e. 


U;, ao? = U;, 1° p 7). (6) 


jon 
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The value of ia, <»» thus obtained (from Uj, .,. 1) may be substituted into equation 
(2) to give the proportional increase of thrust. 


For the stated physical assumptions, this method of calculation is valid and 
exact whatever the values of all the various parameters such as U,/U,, n, p/p; and 
U,/U;. For restricted ranges of these parameters the following argument applies. 


As a preliminary, equation (2) is rewritten as 


Ua, do? 
Ti-Te _ (Waren 1 (7) 
We Ua, oor (N=0) 


(n 1) Ua, 


Substitution of equation (4) then gives 


Ud» cor 
T; T. Uads isn 
Ua, ace (n=0) 


Uas 19 (n=0) 


which can be written more conveniently as 


(n=0) 


The value of (uy. ,/Uas ,) is obtained for general values of n by making the following 
substitutions in equations (5) and (6): — 


U;, ao? + Ua, x? alt U;, (U, + Ua, a) 


p= (5-1): © 


The quantity required is now I’, which is found to satisfy 


r’. U;, 1 Pj Uas (10) 
1+—" (n- | 
[1+ 
U, i, ] 
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If the term represented by ¢, is replaced by a suitable approximation, equation (10) 
may be substituted in both the numerator and the denominator of equation (8) with 
n being put equal to zero for the denominator. This procedure, as compared with 
using, say, equation (3) for the denominator, facilitates the subsequent manipulation. 


For the first approximation it is assumed that ¢, may be neglected as then, if 
BU,/U;(={U,-U,}/U;) is small, the foregoing equations may be solved 
approximately to give the result of equation (1), i.e. 


The only algebraic approximations used in thus reaching equation (1) are those 
subsequent to equation (11). Examination of these approximations shows that 
equation (1) is valid (in the sense of being reasonably accurate) provided that both 


(i) U; mt is small—say, not exceeding 1/10, and 
U, 
(ii) 1} and are not simultaneously large, say 
U, 
.* 1} not greater than 1:5 when n is greater than 3. 


Numerical testing of equation (1) against the algebraically exact calculation, 
using the full equations (2) to (6), confirms these conditions of validity, some 
examples being shown in Table I. The validity conditions would be likely to be 
satisfied in most jet flap applications, not only at take-off but also at cruise, provided 
that the amount of mixing has not been increased above the natural level. 


TABLE I 
COMPARISON OF THE VARIOUS ALGEBRAIC APPROXIMATIONS 
Example p U, ( T,-T, ) 
1-0 10-0 2-0 30 0-208 0187 0-174 
2 3:0 2:0 0:8 1:0 0-033 0-032 0:027 
3 3-0 10-0 2-0 0.162 0-141 0-120 
4 3-0 15-0 15-0 0-600 0°459 0-429 
T,-T, 
Note: ( ; ) represents the result from equation (1) 


T,-T 
( *) represents the results from equation (12) 
2 


T, /fepresents the result from the *‘ exact’’ calculation, using 
. equations (2) to (6) 
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0-20 
CHOROWISE 
STATIC SUCTION PRESSURE 4 
PRESSURE(N. He0.) DISTRIBUTION. 
Ts 
A 
\ ° 
MIXING REGION. 
NOMINALLY SYMMETRICAL AEROFOIL 


AIR BEING ABSORBED (b) 


BY THE JET Cen 
CONDITIONS TUNNEL, OFF 


JET ACTS AS A “SINK” JET NOZZLE DYN. HEAD #30" Hg. 


Fig. 1. 
The sink effect of a jet in “ still ” air. 


The second approximation is obtained by retaining the principle components 
of ¢, in equation (11). This leads to 


Second U 
Approx, 
pp (1 


and Table I shows also the results from this approximation. 


When the validity conditions are not satisfied the full equations must be solved. 
This is straightforward, but no general trends can be seen by inspection. 


The likely significance of the results just derived has been discussed in Ref. 1*. 


3. The Sink Effect of the Mixing 


In Ref. 1 a picture of the flow was based upon the sink effect of the mixing. 
This concept will be elaborated upon and then considered quantitatively. 


*In connection with calculations on the mixing rate, and hence the value of n, the author 
is indebted to C. H. E. Warren for indicating experimental investigations not appreciated 
when considering the value of n in Ref. 1. These investigations were carried out by 
Thompson and Barnes, Squire, Rolls-Royce Ltd., and others, and concerned rates of jet 
mixing when the external stream velocity was non-zero and when the fluid densities were 
unequal. The conclusions adopted in Ref. 1 appear to be consistent with the results of this 
earlier work, 
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TUNNEL SPEEO 100 FT./SEC. 


AEROFOIL SYMMETRICAL, 


JET UNOEFLECTED, 
NOZZLE ALONG TE 


WITHOUT JET WITH COLD JET 
P DYNAMIC HEAD « 
22 IN. 
\ 


4 


static 


(FIG. TO same scace ) 


QUE TO THE MIAING WITH 50% 757 


STREAMLINE 
THE JET ALONG XY 
WHEN NO TET CHORDAL STATION 


THE MIXING ALONG AY CAUSES THE SAME -2-0 (b) EPFECT ON THE AEROFOL 


CONTRACTION CF THE EXTERNAL FLOW AS 


DISTRIBUTION; EVIDENCE OF THE SINK 
WOULD A SINK OF STRENGTH U, (-§w) 


AND OF A DRAG INCREMENT 


Fig. 2 


The sink effect of a jet in a non-zero main stream flow. 


When a jet mixes with surrounding air which would otherwise be stationary the 
jet entrains the air, that is, it draws some of it into itself. To do this the jet must 
create a suction along its boundaries—for example, for the surrounding air to enter 
the jet with velocity v there must be created a suction equal to }pv’, relative to the 
external stagnation pressure at infinity. Now for a pure jet flap having, say, the jet 
undeflected, this suction would extend to the rearward surfaces of the aerofoil but, 
due to weakening with distance, its effect on the leading surfaces would be quite 
small. As a result there would be an increment of drag, the total force experienced 
by the aerofoil being slightly less than the direct thrust from the jet. The difference 
is typically of the order of one per cent. The flow pattern of such a system is 
sketched in Fig. 1(a) and an experimental pressure distribution is shown in Fig. 1()). 


Since the jet is absorbing the external fluid, it is acting, as far as the remainder 
of the field is concerned, in the same way as would a line of aerodynamic sinks. At 
first sight the effect might seem anomalous. The introduction of fluid into the system 
has the effect of extracting fluid from the system; the effect of inserting an aero- 
dynamic source—albeit of a special type—is more readily analysed by replacing 
the source by a distribution of sinks! 


The sink effect of the jet, as this effect is termed, is slightly more complicated 
when the surrounding fluid is moving. Its analysis requires consideration of stream 
tube areas, as indicated in Fig. 2(a). The mixing which occurs along a given portion 
of the jet boundary, say XY in the figure, causes an increase in the width of the jet 
alone, but in general a decrease in the total width of the jet plus the air which enters 
it, this shrinkage being shown in the figure as (- dw). Now, due to this shrinkage, 
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the external flow that has not yet mixed with the jet must be drawn towards it by 
the amount (-4w), and this effect will be referred to as the “contraction” of the 
external flow. However, a two-dimensional line or sheet sink would cause a similar 
flow contraction. Since effectively the contraction represents new boundary 
conditions on the external flow, these conditions being defined uniquely by the 
contraction, and since only the external flow can transmit effects to the aerofoil, 
either from the mixing region or from a sink, it is postulated that there is a sink 
distribution equivalent to the jet mixing and that the criterion of equivalence is that 
the contraction caused by the sinks should be the same as that from the mixing. 
Applying this criterion, and using the fact that a line sink of strength g would cause 
a contraction of q/U,, where U, is the external velocity, the equivalent of mixing 
which causes a contraction (-4w) would be a sink of strength g=U,.(—4w). 
Similarly, for the purpose of a general picture, the distributed sink strength per unit 
area would be U,.(- dw/dx), where w represents the total width of the original jet 
plus a given amount of the adjacent flow. 


The effects on an aerofoil are indicated in Fig. 3. A sink positioned aft of the 
trailing edge of a symmetrical aerofoil in a main stream of velocity U, will cause 


ly 
(a) ZSINKS CAUSE SUCTION PARTICULARLY IN 


THE TRAILING EDGE REGION, HENCE A 


DRAG INCREMENT. 


(b) SEE FIG.3(c) SEE FIG.3@) 
THRUST INCREMENT st DRAG INCREMENT 


PLUS S' FOR SYMMETRY. Wy poTENTIAL FLOW 
MAINSTREAM 
MAINSTREAM VELOCITY VELOCITY DECREASED 
INCREASED AS Ss’ AS ss’ OPPOSE Uo 


IMPLEMENTS Uo ~~ 
~ BOTH EFFECTS Cause 


A THRUST INCREMENT. 
Fig. 3. 
The sinks can induce either a thrust or a drag increment, depending upon position. 
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a drag increment as shown at (a). (Fig. 2(b) shows an experimental pressure 


F 
distribution). A sink placed on the body of an aerofoil, produced there for example d 
as in Fig. 3(b), will cause an increment of thrust. at least for the aerofoil shape m 
shown in Fig. 3(c). 

To utilise the foregoing reasoning quantitatively it is necessary to calculate rn 
firstly, the flow contraction for a given amount of mixing. and secondly, the force - 
which a sink induces on an aerofoil. fl 

Making, as before, the assumption of a square jet profile, the contraction ; 
resulting when a jet of mass flow rate m mixes with n times its own mass flow may | 
be shown to be ; 

U, 
(-dw)=! 
(1 + Ui, 
nu, 
4 
Hence the equivalent sink strength is 
b 
Pi piUj, U, b 


q= (14 
nu, 


It will be noticed that if p,\U,=p,U;, , the sink strength is zero and hence the thrust 
or drag increment is zero. 


A sink of strength g is now considered in potential flow. If the flow is otherwise 
undisturbed and the velocity at infinity is U,, the drag force on the sink is pqU,;: 
this is a standard theorem deducible from momentum arguments. There is also a 
theorem which states that if the flow is that around a body and the sink is positioned 
where the velocity would, but for the sink, have been U,, the force on the sink is 
pqU, parallel to U,. But, from considerations of the overall momentum, the force 
on the whole system must still be pgU,. Hence the force acting on the body must 
be of magnitude pq (U, - U,) (provided that U, is approximately parallel to U,), it 
being a thrust force when U, exceeds U,. This result is readily combined with 
equation (14) to show that the proportional increase of thrust due to the sink effect 
of the jet is 
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For present purposes this gives adequate agreement with equation (1). The main 
discrepancy is the second-order difference resulting from the inclusion of the values 
corresponding to pressure p, in place of those corresponding to pressure py. 


There is a further sink effect occurring when the mixed fluids pass from the 
region of static pressure p, to that of static pressure p,. The change in total flow 
width due to and during the pressure change is different from that for the same 
fluids unmixed. The effect is a second-order one and may be applied as a correction 
to equation (15). Carrying through this correction improves the agreement between 
the results of the various methods, but the algebra of the correction becomes 
cumbersome. Due to the need for this correction, equation (15) is not as accurate 
as equation (1). 


It may be noted that the method just described is related to one sometimes used 
in boundary layer theory. 


4, Mixing Between the Aerofoil Boundary Layers and the Jet 


The mixing between the jet and the boundary layers from the aerofoil may now 
be discussed more fully than was possible in Ref. I. The first five points are taken 
as direct mixing considerations, while the sixth (Section 4.2) is an interpretation on 
boundary layer theory. 


EXTERNAL MAIN STREAM FLOW. 


“SEE * BELOW 


KEY SCALE 


EXPERIMENTAL EDGE OF JET —e— 

[DEPTH OF 0 INDICATES L ! J 
POSSIBLE RANGE OF a) Ye IN. IN. 
INTERPRETATION OF EDGE | 

THEORETICAL BOGE OF JET 


Fig. 4. 


A comparison of the experimental and theoretical jet mixing rates. 
Jet nozzle velocity=48 U,. 


*External boundary layer flow is absorbed by jet in regions AB and A’B’, 
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4.1. Direct MIxInG CONSIDERATIONS 


(i) Equation (21) of Ref. 1, ie. v=(p;/p,)!(U;-U,)/25, shows that the 
entrainment velocity v is proportional to the shear velocity between the jet and the 
surrounding flow. Consequently the rate of mixing will be greater while the 
boundary layer is being absorbed. 


(ii) The oncoming boundary layer will usually be turbulent; after hastening 
the transition from laminar flow at the jet boundary, this would tend to increase 
further the turbulent mixing rate. Moreover, as the boundary layer is likely to be 
much thicker than the jet, its “ mixing length ” will be much greater, but its velocity 
gradients would be smaller. Although the net effect is not clear from purely 
theoretical reasoning, the discrepancy illustrated by Fig. 4 (reproduced from Fig. 2 
of Ref. 1), a discrepancy which exists after allowance has been made for the 
increased shear stress mentioned in point (i), suggests that there is a further increase 
in the mixing rate while the boundary layer is being absorbed. 


(iii) Even for a given value of the entrainment velocity v, the entrainment 
angle tan-' v/u (u being the local velocity in the boundary layer) will increase as 1 
decreases, and will therefore be larger for the boundary layer than for the main 
stream; this would cause increased secondary loss. 


(iv) If the boundary layer has just experienced a pressure rise, its velocity u 
would be still smaller, thus emphasising the increased secondary loss. On the other 
hand, an improvement would be obtained by causing the mixing between the jet and 
the boundary layer to take place in a region of high suction immediately after a 
favourable pressure gradient. This would result in a larger velocity component u 
parallel to the jet, possibly a lower turbulence level, a lower mixing rate, a doubly 
decreased angle of entrainment and_ therefore, probably much _ smaller 
secondary losses. 


(v) The previous four points have considered the effects of the boundary layer 
on only the rate of entrainment, the angle of entrainment, and the secondary loss. 
The boundary layer also affects the primary change in thrust. 


Because the velocity at any position in the boundary layer is less than the main 
stream value, it is more sensitive to changes in pressure. This is shown by 
equation (12) of Ref. 1, i.e. du= - 5p/(pu). As a result the change in the difference of 
velocity (between the jet and the external flow), a change which for the main stream 
was (equation (14) of Ref. 1) 


16) 
pu, 


-(U,-U,) (1 


=0 when pU; 


u pU, 
U, pU 
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It will be seen not only that the change has been magnified but also that there is no 
longer an effective critical condition at which it becomes zero. Thus, returning to 
the energy balance argument of Ref. 1, if the mixing with the boundary layer occurs 
in a region of suction it will cause a primary increment of thrust which is always 
positive, even at cruising speeds. Moreover for an increase in the suction the 
primary gain would increase but, as discussed in point (iv), the secondary loss would 
probably decrease. Consequently a net gain would be expected from the mixing 
with the boundary layer, despite the secondary losses, provided that a local region 
of high suction is practicable. 


42. THE POSSIBILITY OF THE REMOVAL OF FORM DRAG 


The gain of point (v) in Section 4.1 can be re-interpreted on boundary layer 
theory as a reduction, or reversal, of the form drag of the aerofoil. The form drag 
(which is that part of the two-dimensional drag transmitted by the normal pressure 
forces and not by skin friction) may be associated with the boundary layer displacing 
the main flow from the true profile of the aerofoil and thereby preventing the full 
pressure recovery at the trailing edge. 


The rate of growth in the boundary layer of the momentum deficiency (pU ,”6) 
can be expressed by the momentum equation as 


d 
Now the drag of the aerofoil is equal to the momentum deficiency in the wake 


at infinity, thus 


(pU,” wake 


4pU,2c 


(19) 


Integration of equation (18) along both surfaces of the aerofoil therefore gives 


T.E. x 
Cp= 7, dx+ as dx (integrals along both surfaces) 
0 0 


(20) 


The first term in equation (20) is the direct skin friction force (taking the surface 
force as being everywhere parallel to U,), and so the second term must represent the 
form drag :— 


co co 
= |= 
Cp torm= 2c dx ( = d (both surfaces) 
0 0 
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Equation (21) shows that the form drag is produced by thick boundary layers in 
regions of rising pressure, for then 6* will be large and dp/dx will be positive. This 
combination usually occurs towards the trailing edge of an aerofoil. However, any 
device which can manipulate the boundary layer so that the distribution of 6* makes 
the value of the integral zero or negative would be able to eliminate the form drag 
or even to provide a form thrust. Just as for the thrust augmentation considered in 
Ref. 1, this does not imply a creation of power, or even perpetual motion. 


Two examples may be given before examining the effects of injection. The 
first example holds only in supersonic flow, and is taken from Ref. 3. There it is 
shown that, for a body having boundary layer transition at mid-chord, the differential 
behaviour between the turbulent part of the boundary layer and the laminar can 
provide the net gain in thrust. An example in subsonic flow is that of an aerofoil 
with a constant pressure distribution, other than for a “ pulse ” of suction very close 
to the trailing edge. It is arranged for transition to occur instantaneously at the 
peak of the pulse. The value of the displacement thickness 6* would decrease 
sharply at transition and, provided that the magnitude of the pulse were not too 
great, the displacement thickness would be smaller during the pressure rise than 
during the pressure fall; there would thus again be a negative form drag. The 
decrease in the total energy dissipation has been achieved by making the internal 
boundary layer momentum transfer, or mixing, which now occurs largely at 
transition, take place at a high suction, so that the internal velocity differences are 
much smaller. This example is one which might possibly have practical application 
in its own right, but the present question is whether the form drag can be removed 
by high velocity injection. 


To simplify the problem it is assumed that the velocity of the propulsive jet is 
very high compared with that of the main stream. This implies a low propulsive 
efficiency, but it justifies the neglect of both the small negative displacement thickness 
of the jet and also the small changes in the jet velocity that could be caused by the 
pressure variations around the aerofoil. It also shows the source of the economy. 
With this simplification, suppose that the pressure distribution is such that the main 
stream velocity on the upper surface increases continuously to a peak value close to 
the trailing edge and then, after a short length at this level, decreases almost 
discontinuously to the value actually at the trailing edge; further suppose that this 
is made possible by injecting sufficient of the propulsive jet at the beginning of the 
short region of peak velocity just to reduce the total momentum deficiency to zero. 
After the mixing has occurred, the displacement thickness of the boundary layer plus 
jet will become approximately zero and will remain so for the subsequent sharp 
pressure rise. Thus the integral in equation (21) will have a zero value for 6* when 
dp/dx is positive, so that the non-zero values elsewhere will then make the form 
drag negative. The explanation of the mechnical energy balance here again is that 
the high suction decreases the shear velocity in the mixing region and hence the 
mixing losses are reduced. With a lower jet velocity the gain would be less. 


The foregoing arguments all suggest that under cruising conditions it would be 
desirable to arrange for the mixing between the jet and the boundary layer to occur 
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MIXING AND THE JET FLAP 


at a high suction. This would be possible with a shrouded jet flap, and the total 
drag should then be a few per cent. less than that of a conventional aircraft. (In 
contrast the pure jet flap would have a larger drag than the conventional due to the 
mixing with the boundary layer occurring in the high pressure region aft of the 
trailing edge). The arrangement with a region of high suction would produce an 
aerofoil having certain novel features. The section towards the trailing edge would 
probably be thicker and more highly cambered than are conventional sections, while 
the pressure rise would be concentrated into the last few per cent. of the chord. As 
to the performance characteristics, the reduced form drag means that increase of lift 
coefficient would only slowly increase the profile drag, and consequently the 
optimum lift coefficient would be higher than for a conventional aerofoil; the higher 
lift, together with the lower drag previously mentioned, would produce a double 
increase in the maximum lift/drag ratio. However, the increase in optimum lift 
coefficient runs perhaps contrary to one of the basic aims of the jet flap, namely the 
achievement of a large ratio between the lift coefficient at take-off or landing and 
that during cruising flight. 


5. Conclusions 


Of the mixing between the main stream and the jet, one aspect is understood and 
readily calculated. Denoted the “ primary change in thrust,” this aspect may be 
associated with the change in magnitude of the jet and the main stream velocities 


due to the static pressure in the mixing region being different from that at infinity. 
The “ secondary loss ” is not so well understood and cannot yet be calculated. This 
loss is due to the entrained fluid travelling, locally, in a different direction from that 
of the jet with which it is mixing. 


Under cruising conditions the mixing of the jet with the boundary layer seems 
likely to influence both the optimum shape of the wing profile and its performance 
characteristics. The arguments suggest tentatively that the optimum lift coefficient 
and also the maximum lift/drag ratio could be appreciably increased. A detailed 
investigation seems necessary for the development of a satisfactory understanding 
of the flow. 
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Fundamental Principles and Theorems of 
Thermo-Elasticity 


W. S. HEMP, M.A., F.R.Ae.S. 


(Professor of Aircraft Structures, College of Aeronautics) 


SumMMaARY: Kinetic heating associated with supersonic flight presents new 

problems to the aircraft structural engineer. Some of these can be resolved 

by the use of the theorems of thermo-elasticity. This paper outlines the 

foundations of this subject and derives generalisations of the usual variational 
theorems, which are so often used in structural analysis. 


1. Introduction 


Supersonic flight brings, as one of its consequences, significant rises in tempera- 
ture and significant temperature gradients within the structure of the aeroplane. 
The Theory of Aircraft Structures, which has in the past concerned itself mainly 
with isothermal conditions, must therefore broaden its basis and include these novel 
effects of temperature variation. This means that the place previously occupied by 
the Theory of Elasticity must be taken over by the Theory of Thermo-Elasticity. It 
would therefore seem of value, at this time, to enunciate and discuss the fundamental 
principles and theorems of this latter subject. 


NOTATION 


x, y, Z rectangular Cartesian co-ordinates 
T absolute temperature 
X, Y, Z forces per unit volume 


X forces per unit of surface area 


T’ temperature of surrounding medium at body surface 
h heat transfer coefficient 
t time 

T, initial uniform temperature 


u,v, w displacements in the x, y, z directions 


components of strain 


Received August 1955. A simplified version of this paper was given as a Section Lecture before 
the Society on 19th April 1955. 
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THERMO-ELASTICITY 


intrinsic energy per unit volume 
f entropy per unit volume 

0 heat flow into unit volume 

work done on unit volume 


“free energy” per unit volume 


components of stress 


N 

qmQ=—— 


Gibbs function per unit volume 
Young’s modulus 


Poisson’s ratio 


a coefficient of linear expansion 
(1+o0)(1—2c) 

yera- 
jane. Ea 
tinly P= 
a C,, C, arbitrary functions of T (see equations (12) and (15) ) 

It (c)e-0, (c)j-0 heat capacity per unit volume, at zero strain and zero stress 

tal respectively 


l,m, n _ direction cosines of the outward drawn normal 
dV _ element of volume of the body 
dS ___ element of surface area of the body 
Gz, Gy» dz components of heat flux across unit surface in unit time 
K thermal conductivity 


v distance along the outward drawn normal. 


2. Formulation of the Basic Problem of Thermo-Elasticity 


Consider a continuous solid body referred to rectangular axes O(x, y, z). The 
body will be assumed “elastic,” that is to say, the components of strain will be 
taken as single-valued functions of the components of stress and the temperature 
T (equation of state). The body will be subject to a quasi-static loading by forces 
(X, Y, Z) per unit volume and (X,, Y,, Z,) per unit of external surface. The 
temperature 7’ of the surrounding medium at the surface of the body and the 
“heat transfer coefficient” or “external conductivity” h will be assumed known. 
Both the mechanical and thermal loading may vary with the time ¢, but not so 
as rapidly as to excite vibration in the body. The problem to be solved is the calcula- 
sas tion of the deformation, stress and temperature distribution as functions of the time. 


August 1956 185 


terly 


Ww. S. HEMP 


3. Definition of Strain Components 


The displacement (u, v, w) from an initial state of zero stress and uniform 
temperature 7, will be assumed small. The geometry of the body and likewise the 
density will be assumed unaffected by the displacement. The strain components 


Com are then given by* 
ou ov Ow 
Crs Ox Cyy= Ov ’ Oz 
(1) 
awa aw, 


4. Thermo-Dynamical Considerations 


The macroscopic state at a point in the body is determined by the strain com- 
ponents and the temperature. The stress components, intrinsic energy U per unit 
volume and entropy # per unit volume are all functions of these variables. During a 
small change of state, the First Law of Thermodynamics requires that 


. . QQ 


where 5Q is heat supplied per unit volume to the neighbourhood of a point in the 
body and 5W is the work done per unit volume by the surrounding material. 


The Second Law of Thermodynamics requires that 


This relation is valid since the processes are quasi-static and the state of the body is 
sensibly uniform in the neighbourhood of a point. 


If the stress components are denoted by X,,...., Y:,... . it can readily be 
shown that 


Equations (2), (3), (4) givet 
The “ free energy ” F per unit volume is defined byt 
F=U-To . . . . . . 


and is regarded as a function of the strains and the temperature. 


*Ref. 1, §10. 
fRef. 1, §§61, 62 may be consulted for a more complete derivation. 
tRef. 3, Chap. XVI, §3. 
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Equations (5) and (6) yield 5F = - 


oF 
and so oT | 
OF — ol =Y 


The “thermodynamic potential” or “Gibbs function” G per unit volume is 
defined as* 


and is regarded as a function of the stresses and the temperature. From (5) and (8) 


2G _ 
and so or | 


_ 5. Equation of State for an Isotropic Body 


The usual strain-stress relations} for the thermo-elastic case are 


{X,-0(¥,+Z,)} +a(T-T,).... 
2(1 
Cyz= 
where E, o and 2 are all, in general, functions of T. 
Solution of (10) yields 
X,=(A+2p) +A (Cy, + B(T —T,) 
a1) 


where n=E/{2(1+o)}, B=Ez/(1—20). 


Equations (10) or (11) are the “ equations of state ” for the elastic régime of an 
isotropic solid body. 


*Ref. 3, Chap. XVI, §4. 


tRef. 1, §74. 
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6. Formulae for the Functions of State 


Equations (7) and (11) give 
F=} { (A +21) + Cyy + Cre)? + + Con? + Cry? — — — 4€ rx yy) } 
—B(T —T,) (Crr + Cyy + + C, (T) (12) 


where C, is an arbitrary function of T. 


The first of equations (7) gives 


(Cre + Cyy +€:2)— AC, / aT. (13) 


The function C, can be expressed in terms of the heat capacity per unit volume at 
zero strain (c).<». Equations (3) and (13) give 


and so Jeno dT dT (14 


T, 


0 


where it is assumed that F and ¢ vanish for zero strain at T=7T,. Equations (9) 
and (10) give 


(15) 


where C, is an arbitrary function of T. 


Equation (9) then yields 


ale 


+ aT Y,+Z.)+ . (16) 
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THERMO-ELASTICITY 


Introducing the heat capacity per unit volume at zero stress (c);-, gives 


where it is assumed that G vanishes at zero stress and T=T,. 


The formulae (12) and (13) for F and ¢ give, by (6), a formula for U. Another 
expression for U follows from (15), (16), (8) and (11). The consistency of these 
expressions yields the additional relation 


3E2°T 
(c)=0 — (C)e=0= : (18) 


7. Conditions of Equilibrium 


The components of stress must satisfy the usual equations of equilibrium 
throughout the volume and on the surface* : — 


0X, + eX, , 


(19) 


where /, m, n are the direction cosines of the outward drawn normal. 


If the temperature T is given everywhere, the determination of stresses and 
displacements can be effected by the integration of equations (1), (11) and (19) 
Subject to boundary conditions (20) and perhaps given displacements over part of 
the surface. This is however a very difficult problem except in the simplest of 
cases. Special interest is therefore to be found in the following variational theorems, 
which can be used for approximate analysis. 


*Ref. 1, §§47, 54. 
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8. Generalisation of the Theorem of Stationary Potential Energy 


Consider a body subject to given forces and a known temperature distribution. 
Impose a virtual displacement (6u, 5v, 5w) consistent with the kinematic boundary 
conditions. Equation (19) gives 


where dV is an element of volume and the integration is throughout the body. 


aX 
y 


= +X } av=o 


Green’s Theorem gives 


where dS is an element of the body surface. 


Recalling (7) and (20), and remembering that T is not varied, we find 


| | | rav— J [x =0. (21) 


This generalises the usual theorem*. Here, however, the “free energy” F replaces 
the usual strain energy density. 


9. Generalisation of the Castigliano Theorem of Stationary Strain 
Energy 


Consider the same problem as in the previous section, but now impose a 
variation of the stress components (5X,,...,5Y.,....) in such a way that the 
conditions of equilibrium are still satisfied. Equations (1) and (9) give 


ou OG. ] 
ax t OX, dX,+....|dV=0, 


and so 


5 


Ox oy 
+ {| scav-o 


*Ref. 1, “$115, 
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THERMO-ELASTICITY 
Using (19) and (20) in their varied form, we find 
[| Joav+ | @) 


where the surface integral is taken over that part of the surface upon which the 
displacements are prescribed. Equation (22) generalises the usual theorem*. 
Notice, however, that the strain energy density is replaced by the “negative 
thermodynamic potential” (—G). 


10. Heat Conduction 


The determination of the temperature distribution requires a study of heat 
conduction in the body. 


Flux of heat across unit surface in unit time at any point in the body is a vector 


(q:. Gy» Yz). It is related to the temperature by the empirical law 


oT oT 


where K is the “thermal conductivity ” and is a function of T. 


The flux 5Q per unit volume into a small volume element in time 6¢ is given by 


and so, by (23) and (3), 


The boundary condition at the surface is 


+nq.)=h(T’—T), 


Where v is measured hiatiin the outward drawn normal. 


*Ref. 2, equation (64. 10), 
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If the external forces and temperatures are independent of the time, then 
eventually a state of thermal equilibrium is established, whose temperature is 
governed by 


or, if K is sensibly constant, by 


In general the thermal conduction problem is linked through ¢ in equation (24) 
to the state of stress or strain (see equations (13) and (16)). Fortunately in practice 
the terms in stress or strain in the formulae for ¢ and the difference between (c),-; 
and (c),<» in (18) can usually be neglected. Writing these last equal to c, and using 
equations (24), (16) and (17), gives 


= 
K/ ot 0x? dy?’ (28) 


where once more K is assumed constant. This is the usual equation. 


11. Solution of Problems 


(i) The first step to be taken, in the resolution of the problem enunciated in 
Section 2, is the solution of equation (28) subject to the boundary condition (25), for 
the temperature distribution T. The methods for tackling this stage of the problem 
are studied in standard treatises such as Ref. 4, and will not be discussed here. 


(ii) Assuming then that the temperature distribution T is known, it remains 
to analyse the deformation and the stresses. This can, in certain cases, be carried 
out by the direct method of Section 7, but approximate methods based upon the 
variational theorems of Sections 8 and 9 are probably the only feasible methods for 
tackling most aircraft structural problems. There is no difficulty in doing this, 
since it is clear that the usual methods and assumptions will link up just as well 
with the new theorems as with the old. 
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The Theory and Prevention of Undamped 
Aeroplane Nosewheel Shimmy 


D. WILLIAMS, D.Sc., M.I.Mech.E., F.R.Ae.S. 


(Royal Aircraft Establishment) 


SUMMARY: The mathematical theory of nosewheel shimmy is given, with 
particular reference to twin nosewheel assemblies. It is shown that a 
sovereign remedy for shimmy is to make the castor length greater than 
what is here called the “creep distance,” which in practice is found to be 
approximately equal to the tyre radius. Lateral flexibility of the oleo leg is 
disadvantageous but elastic constraint at the pivot is a good feature. The one 
necessitates an increased castor for stability while the other allows a smaller 
castor. It is also shown how, by the use of a compact linkage mechanism, 
the effective castor length can be made independent of the wheel-leg offset 
and can have any desired value. Model experiments that confirm the 
theoretical conclusions are described. 


1. Introduction 


Some time ago the author became interested, as the result of a personal 
experience, in the problem of the snaking of towed caravans, and set out (in a 
private capacity) to find the cause and cure of this phenomenon. Having done that 
successfully he was therefore much interested in the recent crop of “shimmy ” 
trouble which the advent of the twin nosewheel has introduced. 


The problem is somewhat different in the two cases and, of the two, the 
caravan problem is perhaps the more complex. Much of the analysis in the trailer 
problem is, however, directly applicable to the nosewheel problem and there are 
distinct advantages in approaching the nosewheel problem via that of the trailer. 


A main advantage is that with the trailer the pivot point (i.e. the tow-point) is 
far forward of the wheels—the castor length, in other words, is many times greater 
than the radius of the road wheels. This helps to focus attention on the really 
significant characteristics of the tyre, and prevents it being distracted by the host 
of minor but extremely complex forces and couples that flourish within the area of 
tyre-ground contact. When the pivot point is well forward, the inertia moment 
about it is so large that it is quite legitimate to neglect all tyre forces that do not 
make an appreciable contribution towards resisting such moments. On the other 
hand, when the castor length is small compared with the tyre radius, and 
particularly when a single wheel is used, even small forces are not to be neglected, 
because the inertia moments, being practically only that of the wheel about its own 
diameter, are themselves small. 


Originally received October 1954. 
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The legitimacy of relying on an analysis that is based on a fair length of castor 


depends, however, in the last resort on the fact that the final cure here suggested is, 
in effect, an increase of castor length. 


Before starting the analysis of the problem the author wishes to record his 
indebtedness to the work of Temple“ in this country and of Kantrowitz® in 
America for the light they throw on the more recondite aspects of the shimmy 
problem. Mention also should be made of the work of de Carbon in France, 
which incidentally was not seen before the completion of the present work. 


NOTATION 


C,A,, 


= 


t 
l 


x 


x; 


=rifl 


lateral displacement of tyre relative to rim 

forward travel of wheel 

angle made by the wheel rim with its direction of travel 
constant coefficients 

ds/ dt 

time 

horizontal distance between wheel hub and castor axis 


lateral displacement of wheel axle from its central position (or 
undisturbed plane of travel of the wheel) 


lateral displacement of castor axis from its central undisturbed 
position 


radius of wheel (to outside of tyre tread) 
“creep distance ” (as shown in Fig. 5) 
(m+ m,), total mass of system 


first of two masses representing the inertia of the system in 
oscillations about a vertical axis, located at wheel axle, distant h 
aft of c.g. of whole system 


second of two masses above mentioned, located at distance /h, 
forward of c.g. 


v(hh,), radius of gyration of system in plan 
1-(h+h,) 

force per unit lateral deflection of tyre relative to rim 
force per unit lateral deflection x, of the castor axis 


torque per unit angular displacement of the wheels from the 
undisturbed fore-and-aft direction 


m,/m 
F/m 
T /(ml*) 
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A _ root of frequency equation 
X,X,,0 _ initial boundary values of x, x,, y 
Coefficients of frequency equation 


T,,T;,T, test functions 
I moment of inertia of system about castor axis 


q=P/m 
B= 


2. Mathematical Analysis 


Any difficulties in solving problems of dynamic stability are usually connected 
with uncertainties concerning the magnitude and character of the forces associated 
with the various displacement parameters and their derivatives. The main difficulty 
in studying the behaviour of an aeroplane nosewheel is in ascertaining the precise 
character of the action that takes place between the tyre and the ground. 
Flexibility of the wheel-leg in bending and torsion, inertia and gravity forces are all 
easily defined and readily find their appropriate place in the equations of motion, 
but the forces between the tyre and the ground are in a very different case. Apart 
from elastic forces there are, at tyre-ground contact, damping forces and couples that 
are far from being precisely known. 


In the present treatment such damping forces are omitted, in the belief that if 
a cure can be obtained for shimmy trouble by considering the motion of the 
undamped system only, that same cure will be equally effective when damping forces 
are operative. This omission of damping forces is made in the full awareness of 
the important part played by artificial damping in current methods of shimmy 
prevention. Indeed, the whole object of the present paper is, if possible, to deter- 
mine the conditions under which such “ brute force ” methods may be avoided, and 
shimmy prevention achieved rather by reliance on the introduction into the 
undamped system of inherent dynamic stability. 


The first task is to obtain the kinematic relation between the lateral deflection 
of the ground-contact region of the tyre relative to the rim and the forward travel 
of the wheel. Fortunately this has been done and the results placed on a sound 
experimental basis by Kantrowitz™. Those kinematic relations considered relevant 
for the present purpose are presented here and derived from first principles. Only 
one of Kantrowitz’s relations is omitted and that is the one to which he has given 
the name “ kinematic shimmy.” When a tyre is pulled sideways so that it is no 
longer square under the rim, the centre line of the tread where it is in contact with 
the ground is not straight but curved, with the concave side of the curve facing the 
trim. If now the wheel is rolled forward in such a way that, while the mid-point of 
the axle moves forward in a straight line and is constrained against any lateral 
displacement, the axle itself is allowed to rotate about the vertical axis through its 
mid-point, the tyre-ground contact traces out a sinusoidal path. 
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Fig. 1. Fig. 2. 


This happens at dead-slow speed and is independent of inertia forces—hence 
the appellation “kinematic.” The magnitude of the tyre forces involved in this 
kind of motion are, it is thought, small enough to be neglected in comparison with 
the forces considered here. 


2.1. TyRE-GROUND KINEMATIC RELATIONS 


The only tyre displacement that is considered in the present analysis is the 
lateral displacement y of the tyre-ground contact area relative to the rim. In 
Fig. 1, AB is the trace which the vertical mid-plane of the wheel-rim makes with the 
ground, ab is the ground contact of the tyre mid-tread line when undisplaced, and 
a’b’ is the same contact displaced a distance y to the right, being the distance of the 
centre of the contact line a’b’ from the vertical plane AB. 


We are concerned with the manner in which the variation of this distance with 
forward travel (i.e. dy/ds) is affected, first by the angle @ which the rim plane AB 
makes with the direction of motion CD (Fig. 2), second by the first derivative of 
this angle with respect to the forward distance s travelled (d6/ds), and third by the 
value of y itself. These are the only three quantities that are accepted as deter- 
mining the rate of growth of y with forward travel s. 


Since only small values of y are being considered, we may thus write dy/<ds 
in terms of three constant coefficients A,, A, and A. 


dy dé 
Thus ds 7419+ + Asy (1) 


Also, since dy/ds is linear in the three variables 6, d@/ds and y, we may 
conveniently determine the value of the coefficient of each variable under the 
condition that the other two variables are kept at constant zero value. Thus /o 
obtain A, we assume d6/ds and y to be zero and therefore 6 constant. Thus in 
Fig. 3, AB, the plane of the wheel, makes a constant angle @ with CD, the direction 


196 The Aeronautical Quarterly 


of 
unc 
wh 
if 
fur 
P 
the 
anc 

r 
mi 
ax 
thi 
As 
Au 


nce 
this 
ith 


NOSEWHEEL SHIMMY 


| 
| 
| 
| 
b | 
| 
at 


//| pl 


Fig. 3. Fig. 4. 


of travel. The line ab marks the tyre ground contact, and its mid-point P is 
undeflected from the centre plane AB of the wheel. It is convenient to imagine the 
wheel kept stationary and the ground swept back along DC with constant velocity v. 
{f now we assume a small backward displacement ds of the ground and assume 
further that the point P maintains contact (without any skidding) with the ground, 
P comes to P’, and in doing so has moved a distance dy away from the plane AB of 
the wheel. 


It is seen at once that, for small values of 9, 


ds 
and therefore A,=1 in equation (1). 
Moreover, since ds=vdt . (3) 
1 dy 


To obtain A, of equation (1) it is assumed that @ and y are both zero, i.e. the wheel 
mid-axle point is for the moment directly behind, and distant / from, the pivot 
axis O (Fig. 4) in the line of travel, and the point P which marks the contact between 
the centre line of the tyre and the ground is directly under the centre of the rim. 
As the ground moves aft a distance ds, 4 increases from zero to 


dé 
= ds, 


| August 1956 197 


dy 

/ 

A [eae ds 

ds 
the 
In 
the 
ind 
the 

ith 
AB 
of 
the 
er- 
ds 
(1) 
ay 
he 
10 
in 
rly 


D. WILLIAMS 


with the result that the point P on the rim moves laterally a distance P’P, at the 
same time as it moves aft a distance ds from P to P’, so that 


P’P, ao ds. 
ds 


Meantime the point P on the tyre tread, since it is assumed that no slipping takes 
place, must travel straight back in continued contact with the ground to P’. 


After the common aft displacement ds, the tyre and rim are now separated by 
a distance dy equal to P,P’ or 


Thus the coefficient 


(4b) 


It is to be noted that when the total lateral displacement is represented by x, as in 
Fig. 7, 


Oy dx 


To obtain A, we rely on certain experimental work done by Kantrowitz. He 
found, what in point of fact is common experience, that if a tyre is displaced 
laterally from under the rim, and the wheel is then rolled forward without deviating 
from its original plane (i.e. both 6 and d@/ds kept zero) the tyre creeps back under 
the rim before the wheel has rolled through a quarter of a revolution. Moreover, 
the greater the initial displacement the quicker the rate at which it creeps back in 
terms of forward distance travelled. 


Kantrowitz found in fact that the rate of creep is directly proportional to the 


displacement: in other words that 


ds 
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Fig. 5 shows what happens. P is the centre of the tyre-ground contact line 
when the tyre is in its initial undisplaced position, and PP, lies in the central plane 
of the wheel. While the rim is kept fixed the tyre-ground contact point is forcibly 
displaced laterally through a distance y to P, and, by virtue of the load carried by 
the wheel, stays there. On rolling the wheel forward a distance ds along the ground 
(imagined stationary this time) the point P on the rim moves straight ahead to P’, 
but the point P, on the tyre, in moving ahead the same distance, moves also 
laterally to the left to P, in its desire to creep back under the rim. Kantrowitz 
found that the angle made by P,P, with the central line PP, of the rim is such as 
always to intersect that line at the same point H (whatever the value of y) and he 
also found that the distance PH is approximately equal to the radius of the tyre. 
It will be referred to as the “creep distance” or r. It follows at once that 


dy __ 


ds 


and hence that the coefficient 


Again, putting ds=vdf, 
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Pp 


Fig. 6. 


By substituting the values just found for A,, A, and A,, equation (1) can now 
be written in the form 
dy d 


1 


or, with v dt substituted for ds, 


(where again / d6 stands here for the wheel lateral displacement in time 6). 


2.2. EQUATIONS OF MOTION 


The only forces that now remain are easily dealt with, since they are all 
straightforward elastic or inertia forces. 


In the towed trailer problem it was found very convenient, not only in the 
analysis but also in interpreting the results in terms of the physical quantities 
concerned, to regard the swaying or shimmying body (the whole trailer in the 
snaking problem and the whole twin-nosewheel unit in the present problem) as 
consisting of two masses, one situated directly above the axle and the other at a 
point whose location depends on the fore-and-aft distance between the axle and 
the centre of gravity of the unit. 


Fig. 6 makes the matter clear. The wheel, of which A marks the axle, can 
pivot at O about the axis PP, carrying with it whatever masses are attached to it, 
and can shimmy with it, the latter being represented in the figure by the dotted 
rectangle. The centre of gravity of the wheel (or wheels) and the associated masses 
is at B, distant hl forward of the axle A, / being the distance from axle to pivot. 
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Fig. 7. 


Since we are not taking account of any rake of the pivot axis, / defines also the 
castor length. It is convenient to represent the inertia of the system about the pivot 
axis PP by two masses, a mass m at the axle and a mass m, at E, a distance h,] 


forward of the c.g. position B (and «/ aft of O). The masses and their positions are 
related by the appropriate conditions 


m+m,=M=total mass of system, 
mh=m,h,, 
m (hl)? +m, (hl)? =Mk?, 


the moment of inertia of the system about a vertical axis through its c.g. 


The distance of the mass m, aft of the pivot axis is denoted by <«/, so that 
e=1-(h+h,) . ‘ (9) 


2.2.1. The Elastic Forces 
We consider the action of three elastic forces: 
(i) The force F per unit lateral deflection of the tyre relative to the rim. 


(ii) The force P per unit lateral defiection* x, of the pivot axis (relative to 
the fuselage) at the level of the wheel axle. 


(iii) The torque T per unit angular displacement of the wheels from the fore- 
and-aft direction. 


*This may be due to the lateral flexibility of the leg itself or of its connection to the fuselage. 
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Fig. 8. 


By means of (ii) we can represent any lateral flexibility of the wheel leg, and 
by means of (iii) any torsional stiffness of the leg in combination with any torsional 
constraints that may restrict its free rotation. 


2.2.2. The Dynamic System Defined 
Fig. 7 shows the displaced system in plan. 
Since only small angular displacements are being considered, the variable @ is 


discarded and represented by 


where x=lateral displacement of wheel axle. 


The lateral displacement / d@ in equation (7a) is, however, replaced by dx, so 
that 1 d6/dt= x (not x- x,). 


2.2.3. System of Equations 


There are now three equations of motion in the three variables y, x and x,. 
Referring to Fig. 8 and taking moments of all forces about O we have (using dots 
for differentiation with respect to time) 


On putting =p 


=O 


=R 


ms 
m 

F 

m 

mil? 
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this becomes, after dividing through by m7, 
X+p{x+e(x- X,)} ¢+Qy+R(x-x,)=0 
By virtue of the equilibrium of all lateral forces the second equation is 
mk +m, X,)} +Px,+Fy=0 


or, after dividing by m and putting 


P 


X+p +e (X X,)} +qx,+Qy=0 , . (13a) 
The only remaining equation is (7a) which, however, to take account of the extra 
degree of freedom x, of the pivot, must include an extra term x on the right hand 
side. This gives 


or, substituting (x - x,)// for 4, 


Multiplying through by r, and putting 


(16) 


gives nv(x-x,)+rx-ry-vy=0 ? . (16a) 


After collecting terms in the three final equations, equations (11a), (13a), and 
(16a) can be written in the form 


[X (1 + + Rx] + (1 - Rx,]+Qy=0 
[i(1 + + [xp (1 - ©) +9x,]+ Qy=0 @ 
[rx+nvx] - [nvx,]- [ry +vy]=0 
Since these are linear equations with constant coefficients, the three variables 
can be written as 
x=Xe™", x, =X,e"", y=Te™. 
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On substituting these in (17) we obtain 


+ pe*) A? + R] X + [pe (1 R] X, + QV'=0 
[(1 + pe) A?] X + [p (1 - A? +q] X¥,+O0'=0 
[(rA + nv] X+[- nv] X,+(-rA-v)=0 


Eliminating X, X, and I" gives the frequency equation 


[(1+pe*)A?+R], [pe(1-«)A?-R], Q 
[(1 + pe) A?], [e(1-«)A?+q], @Q (18) 
[rA + nv], [ -—nv], [-(rA+v)] 


The expansion of this determinantal equation gives the following quintic in \:— 
Por’ + + + pA? + + (19) 


where, after dividing through by - 1, 


Po=[p 
Pi=[p (1 


P2=[pQ (1 - +g 
P3=[Qpn(1 
Ps=([qQn+qR]v. 


We can now apply Routh’s® well-known method of determining the stability 
of the motion by examining the roots of (19) in terms of the coefficients p. 


The condition for stability is that all the roots of (19) must have negative real 
parts, and this condition is satisfied if the two coefficients p,, p, and the Test 
Functions T,, T, and T, are all positive. 


The Test Functions may be written in determinantal form as 
Pi Do} T;,= Pi Po 0 | 


Ps P2 Ps Pr 
Ps Ps Ds 
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It is seen at once that p, and p, are essentially positive. Also, for the case being 
considered, i.e. where « is very large and negative (see equation (23) et seq.), both 
T, and T, can readily be shown to be always positive. There remains, therefore, 
the simple condition that 7, must be positive. 


2.2.4. Properties of Test Function T, 


It is at once seen, after expanding the determinant T,, that v? is a factor, show- 
ing that there is no such thing as a critical speed and that the motion is either stable 
at all speeds or unstable at all speeds. What prevents an unstable unit from 
shimmying at very low speeds is presumably the smallness of the inertia forces 
compared with the heavy friction damping between tyre and ground, a factor 
omitted in the present analysis. Also the creep distance r (approximately equal to the 
radius of the tyre), except when it appears in the ratio n=r//, appears only as a 
factor r?, so we only need retain n. 


Simplification of Test Function T, 


In expanded form 


T (DiP2 PoPs) (DsPs— PsP2) + Po Ps (22) 


and if the values given for the p’s in (20) are substituted in (22) the expression is 
somewhat cumbrous. However, by confining attention to the particular problem of 
the shimmying nosewheel assembly, for which it can be said that the centre of 
gravity is virtually directly over the axle, i.e. that 


h=0 
and h,= 00 


the coefficients p, and hence the final form of the test function, can be simplified 
greatly. 


Thus m(hl)=m, (h,])=0 ‘ (24) 
so that m, is a small quantity of the second order ‘ ‘ . . (24a) 


and m(hAl)?+m, (hl)? =I : (25) 


where J=moment of inertia of the system about the axis PP’. 


It follows from (24) that m (hl)? =0 and therefore 


m, (h,])? =1=Mk?=mk? . 


where M=m=total mass of system 
k=its radius of gyration about its c.g. in yaw 
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But since h,=(1 - «), 


m, (h,1)? = m, (1 — PF =1=mk?, 
or = mk? 


and (Fy = B (say) 


Also, by (24a), the product of p with < or with any finite quantity is zero. 


The coefficients p of (20) may now be written as 


Do=B 

pPi=B 
P3=q(1+f)+R 
(q+ Q)+Qq 
(nQ+ R). 


Thus T,, as given by (22), becomes finally 


=-q(1—n)(B+n)-BnQ+nR . (28) 


and the stability of the motion depends on whether the right hand side of (28) is 
positive or not. 


The quantities involved are, after dividing throughout by (27), 


q=square of lateral natural frequency of the twin-wheel system on its leg, 
(=P/m) 


Q=square of lateral natural frequency of the twin-wheel system on its tyres 
(tyre treads stationary on the ground and the oscillation excited by lateral 
pushing of rims at their lowest points), (=F/m) 


T T ee T 


mere 


mE ~ m(P+R) td 


=(square of natural frequency of torsional oscillation of twin-wheel system 
about the pivot point) multiplied by (1 + £) 


_(kY radius of gyration of system about its c.g. \? 
B=(\—) = 

l castor length 
n= — fadius of tyre (approximately) 

l castor length 
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In all ordinary nosewheel units the castor length / is much less than the creep 
distance r, so that‘n=r/1> 1. This makes the first two terms on the right hand 
side of (28) negative and the only way to keep T, from becoming negative is there- 
fore by increasing R, i.e. by stiffening up the leg and its steering controls in torsion, 
or, if the wheels are coupled, by stiffening the common axle tube in torsion. 


If, however, the castor length / can be increased so that n becomes < 1, a way 
of making the system stable at once opens out, for all that need be ensured is that 
the first term is greater than the second; the torsional stiffness can then be any 
value down to zero (free castoring with uncoupled wheels). 


The sovereign cure for shimmy is therefore to reduce n=r// well below unity. 
We have little control over r, but the castor length / is amenable to increase. This 
point will presently be discussed in detail. 


Effect of tyre stiffness parameter Q and torsional stiffness parameter R 


It is not surprising to find in (28) that the torsional stiffness parameter R is a 
stabilising factor, but it is odd to find that the tyre stiffness Q is a destabilising 
factor. It looks as if by making Q very large, i.e. by making the tyre excessively 
stiff, instability must supervene. Yet the one and only cause of shimmy is tyre 
flexibility. 


Some light is thrown on the matter by considering what happens when the 
lateral stiffness parameter q at the pivot (i.e. the lateral stiffness of the oleo leg) 
becomes infinite; in other words when the only elastic forces are those due to lateral 
displacement of the tyre relative to the rim and to angular displacement of the 
wheel assembly about the pivot. The second equation in (17) then drops out and 
the criterion for stability is that the test function 


should be positive. 


The essential condition for stability is that n should be <1. Moreover, if 
n> 1, the motion is unstable, however great the value of Q. The explanation is 
that when Q is very large equation (1) no longer represents the true conditions, for 
the tyre is no longer sufficiently flexible to render the assumption of slip-less road 
contact even approximately true. This explains why a very large value for Q is not 
necessarily to be associated with instability in equation (28). 


What is certain is that, according to (28), if the lateral stiffness of the leg is much 
greater (as it normally is) than that of the tyre, and if at the same time n is well 
below unity (0-7 say), the motion is thoroughly stable. This will be seen from the 
typical example given later. 
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What appears to be equally true is that, if m is only just below unity and the 
leg stiffness low, a stiff tyre can cause instability just as a stiff angular constraint at 
the pivot can cure it. Thus there exists the anomalous situation in which the 
danger of shimmy can actually be reduced by letting down the tyre pressure. 


It is to be noted that the extra castor length due to the tyre-road contact being 
a little aft of the vertical through the axle has been neglected throughout. This 
means that there is a small margin of stability in hand. 


It is interesting to note from Temple’s report“ of 1940 that he was fully aware 
of the possibility of shimmy prevention by making the castor length greater than the 
wheel radius and was only deterred from advocating this solution by the current 
view that (to quote from p. 14 of Ref. 1) “ structural considerations prohibit an axle 
offset greater than the wheel radius.” 


Practical values of the parameters 


At this stage it is useful to examine the magnitudes of the various parameters 
concerned, and for this purpose a particular case of a twin-wheel assembly may be 
considered. For this case the following numerical values were found :— 


q= (12:5)?=156, 
Q= (5)’=25, 


R= (13-5)? (1 + 8)=2,200, 
r 22 
n= (7) = 54 


The expression on the right hand side of (28) thus becomes 


156 (- 3) (11+ 4)—(11 x 4-0) 25 + (4 x 2,200)= - 7,000- 1,100+8,800 . (31) 


The value given for R is based on the assumption of coupled wheels, and 
therefore depends on the torsional stiffness of the torque tube connecting the axles 
and on the fore-and-aft stiffness of the tyres. The actual figure quoted is the 
maximum possible for the particular example, being based on a coefficient of friction 
between the tyre and the ground of unity. 


Even so the stability condition is critical and a slight drop in the value of |: 
from the ideal value of unity means the onset of shimmy. 
Uncoupled wheels 


Coupling the wheels is a highly uncertain method of avoiding shimmy and has 
other disadvantages, particularly the bad effect on the steering properties. 
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Consider therefore the effect of uncoupling the wheels and having a free castor, 
| ie. R=0. The last term in (28) consequently disappears and, if n is>1, T, is 
| negative. To make it positive we must clearly have n < 1 by an amount depending 
on the magnitude of the Q term. To see what this entails the typical quantities 
given in (30) may be taken, in which the “ creep distance ” r is assumed to be equal 
' to the tyre radius. The castor distance /, however, is left open, so that (30) becomes 


q= 156, 

{1+ (#)'} 
2 


For convenience we write 


18\? 20.67 n? 
e-(38 n?=0°67 n’, 


so that T, (x constant) =q(1-n)(8+n)-nQ 
= 156 (1 n) (n+ 0-67n?) 0-67 x 25, 


which on dividing through by n leaves a quadratic expression in n. (ie. 
1:28 -0-42 n—n?). 


Thus 7, changes sign when n=0-95, above which value T, becomes negative 
and the motion unstable. A conservative value might be n=0-7, in which case the 
castor length becomes 22/0°7=31 inches. 


The critical value of n depends principally on the creep distance r, which has 
been assumed to be the radius of the tyre. There is some evidence that it is in 
practice less than the radius, although in the model tests described later it was 
1-37 times the radius. Assuming it to be } of the radius would give a critical n of 
0:92 instead of 0-95. A safe value is still given by n=0-7, whence the castor 
length becomes 


3 x 22/0°7=23:°5 inches. 


3. Suggested Practical Means of Increasing the Effective Castor 
Length 


As explained in the analysis, the tangent to the initial path marked out by the 
tyre-ground contact, as the laterally displaced tyre seeks to creep back under the 
rim, cuts the vertical central wheel plane at a constant distance r forward of the 
ground contact. This constant is approximately equal to the tyre radius, but can 
only be found by experiment (see Section 4). 
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Since the main condition for stability is that the castor distance should be 
longer than r, the value of r is vitally important. It can be found experimentally, 
by marking or colouring the tread centre line and examining the trace left on the 
ground by its path as the laterally displaced tyre creeps back under the rim, the 
latter meanwhile being constrained to move forward in its own vertical plane, as 
described by Kantrowitz. 


The safest plan, however, is to design the wheel assembly to have / comfortably 
longer than r, and this will be ensured if it is made at least equal to the tyre 
diameter. 


3.1. PRINCIPLE OF SCHEME 


To do this by offsetting the wheel from the oleo leg is not very feasible. It is 
therefore suggested that the way to overcome the difficulty is to incorporate the 
castor inside the fuselage as shown in Fig. 8 on p. 202. The leg is represented by AB 
and the under surface of the fuselage by ab. The vertical reaction from the ground is 
taken either at A or A’ or shared between A and A’, and the function of the links 
OA and O’A’ is to ensure that the whole wheel assembly rotates bodily about the 
axis OO’, a distance L, of about twice the tyre radius, ahead of the leg. 


Admittedly this idea complicates the method of attachment of the leg to the 
fuselage and some extra weight cannot be avoided. To counter this, however, 
something is gained by doing away with the external castor, i.e. the offset between 
leg and wheel axle, since the bending moment in the leg introduced by this offset 
is thus eliminated. 


3.2. PRACTICAL APPLICATION OF PRINCIPLE OF SCHEME 


It will in some cases be impracticable, and in all cases inconvenient, to house 
the castor arms OA, O’A’ which, having a length nearly equal to the diameter of 
the wheel, will be something like a yard long. This difficulty can be overcome by 
using a compact linkage mechanism which, while itself extending no more than a 
foot or so forward of the leg, can be so designed as to locate the effective pivot axis 
any desired distance ahead of the leg. 
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UPPER LINKAGE 
£ B 


LOWER LINKAGE 


Fig. 10. 


| 32.1. Linkage Mechanism* 


The principle of this linkage mechanism can be understood at once from the 
; plan view shown in Fig. 9. 


A link AB pivoted to a fixed point A is pin-pointed at B to a second link BE, 
sliding in the end E of the link DE and fastened integrally to the oleo leg at B so 
_ that any angular displacement of the link BE is transmitted to the leg at B. 


A third link DE, pin-jointed to the link AB at C, has one end D sliding in the 
fixed straight groove ab and its other end E sliding along the link BE. 


By suitably relating the ratios BC/CA and EC/CD the line of the link BE 
produced can be made to cut the base DA produced at any desired point O. This 
point is then the effective pivot point, and OB is the effective castor length. 


Although geometrically easy, it pays not to have the point O unduly far 
forward, because the farther forward it is, the smaller the angular displacement of 
the link BE (and hence of the leg) for a given size of linkage (i.e. length of links). 
If the pivot O is located some 36 inches ahead of A, an angle (BOA) of about 8° 
can be reached without more than 6 inches lateral displacement of the leg. 


Large-angle displacement of leg 


The 8° amplitude of leg movement just referred to covers all high speed 
travel of the nosewheel, so that for such speed the linkage is fully operative in 


*The linkage mechanism described here serves only to explain the principle involved. A 
Practical mechanism embodying this principle, but differing radically in design, has since 
been developed. 
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locating the pivot forward. It is, however, necessary also that sharp turns be 
possible at very low speeds when the danger of shimmy is negligible. To cater for 
this condition the linkage is so designed that, once the limiting value of = (about 8°) 
is reached in a cornering manoeuvre of the aircraft, the effective pivot point is then 
automatically translated from O to A (or even to B itself, with a little castor at the 
wheel), with an effective castor thereafter of length AB. This translation is 
achieved by ending the slot ab at D,, where a branch slot FF’ opens out in the form 
of an arc of a circle with centre A. Any rotation of the nosewheel beyond the 
limiting value of z occurs about A as centre and is accompanied by a sliding of the 
end D of ED along the circular arc D,F. 


3.2.2. Design Problem 


Fig. 10 shows in elevation two sets of the linkage mechanism AB and A’B’, the 
leg BB’C being guided by the two arms AB and A’B’. These are constrained to 
rotate together by being carried on the same vertical tubular sleeve which can turn 
about the axis aa. 


The main design problem is to arrange that, while being displaced laterally by 
the rotation of the arms, the vertical ground load is adequately resisted at B or 
B’ or at both. Possibly the lower support B’ may be arranged to travel on a roller 
track. Whatever method is adopted has to be compatible with the ancillary 
retracting mechanism of the wheel assembly. 


This, however, is not the place to discuss design details which must in any 
event be dependent on the individual layout of the wheel assembly. 


4. Experimental Work 


To check the theoretical conclusions described, a small model was made in 
which two model-aeroplane wheels were used to make up a twin-wheel assembly. 
Means were arranged for : — 


(i) Keying the twin wheels to rotate as one or allowing independent rotation. 
(ii) Changing the tyre pressures. 
(iii) Changing the lateral flexibility of the leg. 
(iv) Changing the elastic angular constraint of the system about the pivot 
point. 
Fig. 11 shows a side elevation sketch of the apparatus. 


The pivot bush B,B, is connected to the fixed vertical support AA through 
a parallelogram linkage which is spring-tongued at its corners t. The pivot spindle 
is thus flexibly supported in the line of its axis while motion about any axis normal 
to this line is prevented. The arms C, and C, rotate integrally with the spindle. 
and carry the leg PQS to which the twin-wheel assembly at S is attached. The leg 
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Fig. 11. 


is pinned at P and passes through a slot in the arm C, so that the fore-and-aft 
position of the twin wheel can be set at any desired value of castor length / by 
slackening the butterfly nut at Q, sliding the leg into its new position, and screwing 
up again. 


Forward travel is simulated by resting the wheels on the periphery of a drum 
D which is rotated by hand in the direction of the arrow. 


4.1. EXPERIMENTAL RESULTS 


As far as could be observed, the shimmy behaviour, as the theory indicates. 
was independent of the speed of rotation of the drum. 


First experiment 


In the first experiment the leg was made rigid against lateral flexure and the 
whole assembly was free to rotate about the axis B,B,. The behaviour was found 
for slack tyres and for tyres moderately tight. Thus, in terms of the criterion given 
by (28), 


ie. g(l1-n)(B+n)- BnQ+nR, 


which must be positive for stability, g is large, Q is small or moderate according as 
the tyres are slack or under moderate pressure, and R is zero. 
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TABLE I 
Slack tyres Moderate tyre pressure 
Up to 1-4r, Spontaneous shimmy Spontaneous shimmy 
Oscillation dies Spontaneous shimmy 


down in 3 or 4 cycles 
after finger flick 


j 
| 


1-72r, Stable. Oscillation dies Stable but dying down only 
down immediately slowly 
1:86r, Thoroughly stable Stable but taking few cycles 


to die down 


Table I shows how the behaviour varied with different castor lengths, expressed 
in terms of the actual tyre radius r,. 


Two facts emerge. The first is that the creep distance r, which in this instance 
must be close to the critical castor length, is more like 1-5r, than r,. The second 
is that the shimmy is not sensitive to tyre pressure when the leg is rigid. 


The second fact was to be expected since, in the criterion, the second term is 
small compared with the first for all values of n except very close to unity. 


The fact that the creep distance, judging by the critical castor length, appeared 
to be nearer 1:5r, than the tyre radius r,, raised the question of what the actual 
creep distance was. A simple experiment was therefore carried out on the lines 
suggested in Section 3. 
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The leg of the twin wheel was gripped firmly in a horizontal position in the 
jaws VV of a vice as in Fig. 12, and a vertical smooth surface A was fixed adjacent 
to it. Another vertical smooth surface BB,, to which a sheet of squared paper was 
stuck, was held gently but firmly against the twin tyres so as to leave a small gap 
between the edge B, and the surface A. The face BB, was now shifted laterally in 
its own plane (still maintaining non-slipping contact with the tyres) to make contact 
with the surface A. The tyre contacts with the face BB, were thus forced out of 
the plane of the wheel rims by an amount equal to the width of the gap. With the 
vertical edge B, in continuous contact with the face A, the face BB, was then pulled 
upwards, the tyres meanwhile rolling and the tyre contacts gradually creeping back 
under their respective rims, leaving visual traces of the paths followed by the tyre 
contacts. 


This simple experiment was repeated a few times with practically identical 
results each time. It was found that the central tyre ridge moved laterally from 
0:15 to 0-05 inches towards its central position under the rim in tracing out a for- 
ward distance of 1:5 inches, i.e. 1-5r, (since the tyre radius r,=1 inch). 


Representing the curve by the expression 


where s is the forward travel as in (5), we find the value of p by writing down 


— ps 


Ae 
Ae"? 
84) 

Y2 


2 


Substituting 1:5 r, for (s, —s,) and 


0-15 v1 
0-05 =3 for Yo 
it i : ; 
is seen that log. 3 = 
0 0 


Thus =-py=- (x) 


and, by comparing this with (6), it is seen that the creep distance for the small fat 
tyres of this model is 1:37 times the tyre radius. This now fits in well with the 
critical castor length of about 1-5 r, to 1:6 r, found, since the presence of the second 
term in the criterion expression of (28) requires that the critical castor length must 
be somewhat greater than the creep distance. 
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TABLE Il 
FLEXIBLE LEG, FREE PIVOT 
Remarks 
Castor length 
Slack tyres Moderate tyre pressure 
Zero to 2-7r, Unstable Unstable 
2°8r, Steady mild oscillation of | Spontaneous energetic 
small amplitude shimmy 
3-0r, Stable, unless disturbed eed 
3°4r, Stable Stable 


Second experiment—Effect of lateral flexibility of the leg 


The second experiment was made under the same conditions as the first, except 
that the leg was given a good measure of lateral flexibility, the force at the wheel 
hub per unit lateral displacement being comparable to that required for tyre 
displacement with the tyres in the slack condition, and much smaller than that for 
the tyres in the blown-up condition. Table II gives the results. 


It is seen here that the introduction of lateral flexibility in the leg has made the 
system much more ready to shimmy and that, whereas a castor length of about 
twice the wheel radius ensured stability in the rigid leg condition, a castor of more 
than three times the radius is necessary with the flexible leg. 


If the criterion, i.e. the expression (since R is zero in (28)) 


q(1-n)(8+n)—8nQ 


is now examined, it is seen that this result was to be expected. For, to keep the 
expression positive with the reduced value of g clearly requires a reduction in n. 
As n is reduced the negative second term quickly drops in value whereas the first 
changes little. 


The fact that a greater castor length (i.e. a smaller value of n) is necessary for 
stability with the harder tyre (i.e. with a higher value of Q) is also in agreement with 
theory. 


Effect of elastic constraint of pivot angle 


Although the provision in actual design practice of a stiff elastic constraint at 
the pivot does not seem to be a very feasible idea, it is nevertheless instructive to 
see the effect of such a constraint, if only for the way in which the theory is 
confirmed by experiment. 


To introduce the necessary constraint, the top arm C, (Fig. 11) was attached 
to lateral springs running perpendicular to the plane of the paper and anchored to 
two rigid arms fastened to the fixture E. 
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Two sets of experiments were made, the first with a comparatively rigid leg and 
—___ [| the second with the very flexible leg used in the Second Experiment already 
| described. 


Comparison of the results with theory is facilitated if the criterion expression 


q(1-n)(6+n)—BnQ+nR 
is written in terms of masses and lengths rather than in terms of natural frequencies. 


ure 


Thus, since q= ,Q= R= B= 


and n=r/I, it is found, on substitution, that the criterion takes the form 


except 


r tyre 


at for 
For a given leg stiffness P, a given tyre stiffness F, a given angular stiffness T, 


a given radius of gyration k of the wheel assembly about its c.g. and a given creep 
de the | distance r, expression (33) can only vary with the castor length /. 


about 
more @ Lhird experiment—Rigid leg with elastic angular constraint 


The leg in this experiment was again made rigid and a spring angular constraint 
was introduced at the pivot. This means that, in (33), P is large and T has now a 
finite value instead of the zero value corresponding to the free pivot of the previous 
experiments. As before F was given two values, one corresponding to low tyre 
pressure and one to moderate pressure. 


the 
ge A glance at (33) shows that the expression would be expected to become 
; al positive at a lower value of castor length / than that for T=0. Table III, which 
gives the results for moderate tyre pressure (the lower pressure required less castor 
to overcome instability under all conditions) shows that experiment confirms this. 
ry for 
t with TABLE Ill 
Remarks 
Castor length 
Free pivot With spring constraint 
Otol7 Spontaneous shimmy 
int at 2°1 ry (or >) Completely stable 
ive 0 to 11 4, Spontaneous shimmy 
ry is 
14 1, — Shimmy on slight 
disturbance 
ached 
ed to 18 ry Completely stable 
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TABLE IV 


Spring constraint 
on pivot 


Castor | 
Castor length Free pivot Castor length 


0 to Spontaneous 0°5 Stable 
2°8 ry shimmy 09 r, Slight steady 
vibration 


1:4 r, to Spontaneous shimmy, 
most violent at 2:Ir, 


Stability 2'8 ry Small steady 
increasing oscillation 


Stable 3°2 Stable 


In the Fourth Experiment the spring constraint at the pivot was retained but 
the leg was again made flexible. The results (for moderate tyre pressure) are given 
in Table IV. 


The interesting feature of Table IV is the way in which, with the pivot 
elastically constrained, there is stability at very low and at very high castors but 
instability in the middle range of castor length extending from /=1-4r, to /=2:3r, F 
To explain this behaviour notice that, in the criterion (33), if 


(P+F)<T 


by only a small amount, the criterion must become positive for some sufficiently 
small value of the castor length /. It is also clear that stability must be reached at 
the other end of the scale for a sufficiently large value of /. Whereas stability can 
always be ensured by the latter course, the existence of a stable region at small 
values of //r depends on the value of T in relation to P, F and k’. 


A further test, in which the pivot spring stiffness was halved, did show that [ 
instability persisted right down to low values of the castor /. At zero value of 7, — 
as Tables IV and II show, stability is only to be attained by a comparatively high f 
value of 


Fifth experiment—Coupled wheels 


With the two wheels of the twin-wheel unit coupled together to rotate as one, § 
no shimmy was found under any conditions. This is the only experiment in which ; 
damping was deliberately introduced, and it shows that in this particular model the 
amount of damping was adequate to prevent shimmy. Experience with full-scale 
work indicates however that shimmy persists in spite of wheel coupling. 
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5. Conclusions 


The model tests described are admittedly more of a qualitative than a quantita- 
tive character, but, so far as they go, they do confirm the main conclusions reached 
by the theoretical approach. These may here be conveniently summarised as 
follows : — 


(i) Adequate castor length. An adequate length of castor eliminates the 
possibility of shimmy. 


(ii) Factors governing critical castor length. The precise value of the 
critical castor length depends on the “creep distance” and on the 
relative values of tyre lateral stiffness relative to the rim, the oleo leg 
lateral stiffness and the angular stiffness (if any) of the system about the 
pivot axis. 


Free pivot with stiff leg. With a free pivot, and a leg stiffness (measured 
at the wheel axle) several times greater than the tyre stiffness (measured 
at the wheel axle relative to the ground), the critical castor length is very 
nearly equal to the creep distance—itself normally between 1-0 and 1-5 
times the wheel radius at the tyre tread. 


Free pivot with flexible leg. With a free pivot, and a leg stiffness 
comparable with, or lower than, the tyre stiffness, the critical castor 
length is much increased and can well be greater than twice the creep 


distance. Lateral flexibility of leg under free pivot conditions is there- 
fore favourable to shimmy. 


Tyre stiffness. The effect of increasing tyre stiffness (or tyre pressure) 
is to increase the tendency to shimmy and to necessitate a greater castor 
to stop it, but the system is not very sensitive to the amount of variation 
of stiffness normally available. 


Elastic moment at pivot. An elastic angular constraint at the pivot is 
invariably a good feature and, if a strong elastic constraint is applied 
there, the increased castor length necessary to offset the low leg stiffness 
mentioned in (iv) is appreciably reduced. 


Strong elastic moment with flexible leg. An interesting corollary of (iv) 
is that, although leg flexibility is a bad thing normally, it could be a 
valuable feature if only it were possible to combine it with a strong 
elastic constraint at the pivot. For, under these conditions, the system 
becomes stable at very small castor lengths. Unfortunately, however, 
there seems to be no obvious way of exploiting this feature. 


Coupling the twin wheels. The effect of coupling the two individual 
wheels of the twin-wheel unit so as to rotate as one is to introduce an 
amount of damping that may inhibit shimmy entirely. In the experi- 
ments made on the model discussed here such coupling did inhibit 
shimmy under all the conditions tried. In any particular full-scale 
example the efficacy of wheel coupling would no doubt depend on the 
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degree of torsional rigidity of the common axle and upon the magnitude 
of the inertia forces compared with the tyre-ground damping forces. It 
is, however, to be emphasised that no attempt was made to include 
such damping forces in the present analysis and the experiment of 
coupling the wheels was only included as a matter of interest, since it is 
freedom from shimmy with uncoupled wheels that is of primary interest. 


Mechanical device for increasing castor length. It is possible, as 
described in Section 3, to shift the pivot point forward artificially so as 
to provide the length of castor necessary for stability. The linkage, or 
gear device, required for this need only cater for the small castor angles 
associated with forward travel at speeds above the comparatively low 
manoeuvring speeds that require large castor angles. For such large 
angles the pivot is arranged to automatically shift back close to the 
vertical line through the tyre-ground contact point. 


(x) Special small wheels and/or special tyres. Since what is required for 
immunity from shimmy is a large value of the ratio //r=(castor 
length)/(creep distance), such immunity can be attained by reducing the 
creep distance just as well as (and in fact to much better practical 
advantage than) by increasing the castor length. This reduction would 
follow a reduced nosewheel radius or the use of a tyre with a creep 
distance as low as possible in relation to the tyre radius. Exploration of 
the latter possibility may well lead, in the author’s opinion, to 
valuable results. 
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Vibration of Isosceles Triangular Plates having 
the Base Clamped and Other Edges 
Simply-Supported 
HUGH L. COX and BERTRAM KLEIN 


(T. Y. Lin and Associates, Los Angeles, California) 


SuMMaRY: Approximate solutions are given for the lowest natural frequency 

of flexural vibration of isosceles triangular plates which have the equal edges 

simply-supported and the base clamped. The results are presented graphically 
to allow rapid determination of a desired fundamental frequency. 


1. Method of Solution 


The approximate solutions are obtained by using the method of 
collocation“: *», Let x and y be co-ordinates in the middle surface of the uniform 
elastic plate shown in the lower part of Fig. 1. If the static loads on the plate are 
replaced with equivalent inertia forces when the plate is vibrating freely, the 


governing differential equation of the plate may be written as 


| where D is the plate stiffness, w represents the plate deflections, positive downward, 


q is the weight of the plate per length squared, g is gravitational acceleration, , is a 
natural frequency and V? =0?/0x? + 0?/dy’. 


The boundary conditions are 


Ow 
(=) 3) 


=0 ‘ (4) 


z=+(a/h) 


where n is the normal direction to the lines x='+(a/h) y, and ¢ is the tangential 
direction of any line along a rectilinear edge. The term 0*w/@?* is equal to zero. 
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Vibration coefficients for triangular plates having the 
base clamped and equal sides simply-supported. 


The deflection function w= cos( x). 


ny *Y sin 272 


where F,= y? sin » sin » (y h)’], 


and z=generalised coefficient, h=altitude of triangle, a=half of base of triangle, 


satisfies equations (2) and (3) exactly and the condition 


On? 


z=+(a/h) 

from equation (4) at some point in the region h/2<y<3h/4. Selected points 
for collocation will be on the line x=0; therefore, no serious error should result 
from not satisfying completely the boundary conditions. The first term of 
equation (5) appears physically to have the shape that one might expect for a 
triangular plate vibrating in its fundamental mode. 


It will be necessary to differentiate the deflection function given in equation (5) 
and substitute the proper derivatives into equation (1). Let 
F (y)=2,F, + + 25F;. 


Equation (5) may be written w=F (y) cos > *) 
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VIBRATION OF PLATES 


Values of w and derivatives of w at x=0 are 
(w):-0=F (y) 


(2) 


_ @F(y) 

OF"? dy* 


If derivatives of F (y) with respect to y are taken and are substituted into equations 
(8), and if equations (8) are then substituted into equation (1), there results 


Az,+Ba,+C2z,=0, at y=h/2 
Ja,+Ea,+F2,=0, at y=2h/3 
at y=3h/4 
4 2 
42-0523+97-4100(") - 0-258 
a a 
B= 1736-351 
4 2 
c=1-52202 (*) - 6+4-9348 0-0156258 
a a 


-0-333338 


E 10-2736(")' 1130-014 231-347 


vints F=0-67645 + 24+ 3-28988 Cy 00219478 
esult 
4744-15 - 0281258 
Or a a 


H= ~1-65301(")' 2879-837 -172:190(") +0:397758 


I= 038051 (/ + 46: 500+ 1:23370(" 4] 0:0197758 
(6) 


(7) and where 
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Values of 8 for various ratios of h/a may be determined from the condition 


ABC 
J EF\=0 . ‘ ‘ » 


Fig. 1 shows in graphical form, on semi-log paper, the relationship between 
the vibration coefficient y and h/a, where y=}, for a practical range of h/a values. 
For the particular case of a right-angled isosceles triangle, from equation (12), 
y=34-7. The fourth mode of flexural vibration of a clamped square plate 
corresponds to this case. An approximate answer for the fourth mode has been 
found by Young to be y=32-9. There is about five per cent. difference between 
the two solutions. The collocation solution is believed in this case to give answers 
that are slightly too high; however the curve in Fig. 1 affords a rapid approximate 
answer that should be reasonable for engineering purposes. Certainly the results 
are sufficiently accurate for aircraft panel studies, since the boundary conditions 
seldom are known exactly. 


If the plate shown in Fig. 1 rests on a uniform elastic foundation, then the 
fundamental frequency may be written® as 


where the symbol K represents the elastic foundation modulus per length cubed. 
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On a Class of Laminar Viscous Flows Within 
One or Two Bounding Cones 


A. J. A. MORGAN, Ph.D. 


(Aeronautical Engineering Research Inc., Pasadena) 


Summary: The class of incompressible axially-symmetric laminar viscous 
flows considered are exact invariant (similarity) solutions of the Navier- 
Stokes equations. The functional forms of the velocity components are 
deduced by group-theoretic arguments. The system of governing partial 
differential equations is reduced to a system of ordinary differential equations 
and these in turn are reduced to a single first order non-linear ordinary 
differential equation for the dimensionless tangential velocity component. 
The general solution of this equation is obtained in terms of hypergeometric 
functions. 


The classes of laminar viscous flows bounded by a right conical surface 
or by right inner and outer conical surfaces with a common apex obtainable 
from the previously mentioned (similar) solutions of the Navier-Stokes 
equation are studied. Suction/injection velocities (inversely proportional to 
the distance from the apex) normal/parallel to the conical boundaries are 
admissible for the classes of flows considered. If the flow is bounded by 
impermeable conical surfaces, then it is shown that these solutions of the 
Navier-Stokes equation imply that it must be identically quiescent. 


When the bounding cones are not impermeable, closed solutions for the 
single-cone case are found in terms of elementary functions. A numerical 
example of such a flow, with fluid injected normal to the conical boundary, 
is given. Such a simplification of the hypergeometric functions entering into 

the problem, however, cannot be achieved in the two-cone case. 


1. Introduction 


The class of incompressible axially-symmetric laminar viscous flows considered 
can be characterised as those in which the component velocities are inversely pro- 
portional to the distance from the origin. These flows are of interest since they 
constitute one case in which the Navier-Stokes equations can be integrated exactly. 
This latter fact has been previously and independently recognised by Yatseyev'’) and 
Squire), Yatseyev initiates his analysis by assuming that the velocity components 
are of this form; Squire, by assuming the equivalent special form for the stream 
function. In the present work this form for the velocity components is deduced by 
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showing that it is the only one which will yield invariant (similarity) solutions*) of 
the applicable continuity and Navier-Stokes equations under the one-parameter 
continuous group of transformations which is the generalisation of the group of 
uniform expansions. This aspect of the problem, together with a discussion of the 
boundary conditions applicable at the surfaces of the cones for this class of flows, is 
considered in. Section 2. 


Since the class of flows under consideration is described by invariant solutions 
of the applicable continuity and Navier-Stokes equations, they can be transformed 
into a system of ordinary differential equations which, in turn, can be reduced to a 
single first order non-linear differential equation, containing three arbitrary constants 
of integration, for the dimensionless tangential velocity component. Yatseyev'') has 
obtained the general solution for this latter equation. For the sake of completeness 
a brief account of a somewhat different method of deducing this solution is given in 
Section 3, since some new relations between the dimensionless velocity components 
are obtained thereby. These relations are of importance in the present study of 
viscous flows within one or two bounding cones. 


Several special solutions of this first order ordinary differential equation have 
been obtained by making restrictive assumptions on the magnitude, or inter-relations 
between, the three arbitrary constants of integration which it contains. In this 
manner Yatseyev'') and Squire, in somewhat greater detail, have discussed the 
laminar flow produced by the continuous application of a force at a point in a 
viscous fluid of infinite expanse. This latter fundamental description of the flow is 
due to Squire, who descriptively refers to the flow as that due to a round laminar 
jet. Yatseyev” refers to this same flow as that due to a stream flowing out of a thin 
pipe into an infinite region filled with the same fluid. According to Yatseyev, the 
solution to this particular problem was originally found by Landau™. Under 
another set of assumptions on these constants of integration, Squire“ has discussed 
a viscous flow which he classifies as that due to a “ jet emerging from a hole in a 
plane wall.” In view of the present work a different interpretation for this flow will 
be given herein (Section 7). Under still another set of these assumptions, Yatseyev"’ 
has discussed the flow from a half-line source or sink into an infinite region filled 
with the same viscous fluid. These special cases, to the author’s knowledge. 
constitute all of those which have so far been covered in the literature. 


Slezkin“ did not (as otherwise implied in Ref. 1) obtain or investigate the exact 
solution of the Navier-Stokes equations considered by Yatseyev’”, Squire’*’, and the 
present author. 


The possibility of studying those solutions of the first order ordinary differential 
equation which describe laminar viscous flows within one or two bounding cones 
does not seem to have been previously recognised. It appears, in fact, that this is 
the naturally admissible interpretation for such flows, since it gives a means of cor- 
relating the four arbitrary constants of integration contained in the general solution 


*Descriptively speaking, those solutions of partial differential equations whose functional form 
is such that the number of independent variables in these equations can be reduced by one. 
See Ref. 3 for precise definitions, 
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* of of this equation with the pairs of suction/injection velocities which can be assigned 
neter at each of the two conical boundaries. Such an investigation is pursued in the 
ip of present work. The admissible suction/injection velocities are inversely proportional 
f the to the distance from the apex and may be parallel/normal to the conical boundaries. 
NS, IS The mass flow rates for flows of this class existing within one or two bounding cones 
are discussed in Section 4 and the detailed description of such flows is studied in 

tions Sections 5 and 6. 
&rmed 

toa NOTATION 
Vhs a _ parameter in a continuous one-parameter group of transformations 
ness b parameter occurring in solutions for 7 (€), see equations (31) 
sents C constant of integration occurring in equations (18) or (34) 
y of C,,C., constants of integration occurring in equations (31) 
F(z, 8; y;&) hypergeometric function with parameters z, 8, y, and argument € 
ions K,L,M _ constants of integration in equation (13) 

this 

the M _ mass flow rate across that portion of a spherical surface, centred at 
a the origin, between and intersected by the two bounding cones 
w is m,q,r,S exponents associated with parameter a, see equations (4) 
thin P (6) = - pR*/v?, the dimensionless pressure at a point in the fluid, a 

function of 6 only 

the 
nder §— Pp pressure at a point in the fluid, a function of R and @ due to 
ssed imposed axial symmetry 
P.q exponents in the transformation connecting ¢(€) and »(£), see 
will equations (17) 


led — Pq; i, j=1, 2; particular values of p and q given by equations (19) 
and (20) 


dge, 


spherical co-ordinates 


oo t(€) transformation for V (4), defined by the first of equations (17) 

the 

° U (@) = uR/», the radial Reynolds velocity component, a dimensionless 
function of 6 only 

tial 

U;,U, values of U (6) at 6=A and A respectively 


nes 
sis § V (#) = vR/», the tangential Reynolds velocity component, a dimensionless 
function of only 


ViiV, values of V (6) at and A respectively 

u,v,w radial, tangential (i.e. in 6@-direction) and azimuthal (i.e. in 
¢-direction) fluid-velocity components 

a,8,y parameters of hypergeometric function F (z, 8; y; €) 
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i, j=1, 2; parameters of hypergeometric function occurring in 
solutions 1, (€), given in terms of p;, q;, and K by equations (21) 
and (22) 


2p;, third parameter of hypergeometric function occurring in 
solutions (€) 


transformation for ¢ (€) defined by the second of equations (17) 


a second transformation for V (6), defined by the first of 


equations (29) 


solutions of the transformed equation (14) associated with the pair 
(pi, qi) 


semi-vertex angle of outer bounding cone 
semi-vertex angle of inner bounding cone 
kinematic viscosity of fluid, a constant 
cos? (4/2) 

values of € at 9=A and A respectively 


density of fluid, a constant. 


2. Formulation of the Problem 


The spherical co-ordinate system (R, 6, ¢), the corresponding velocity com- 
ponents (u, v) and the geometrical aspects of the problem to be solved are illustrated 


in Fig. 1. 


In this co-ordinate system the Navier-Stokes equations, on suppressing 


the third velocity component w, can be written as 


OR 


“OR Roe 


uv 
~ aR? * 


Oru 2 ou 1 


2 cotddu 2u 
OR* RoR  R*0e 


R’? 06 


1 op 


2 dv cotday 

ROR R* 00? 06 
R? 06 ~=R* sin’ @ 


1 op 


and the continuity equation becomes 


where p is the static pressure, and p and v denote the density and kinematic viscosity 


of the fluid respectively. Also, due to the imposed axial symmetry, 


u, v, p are 


functions of R and ¢ only. 
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Fig. 1. 
The two bounding cones and co-ordinate system. 


Only that class of flows given by the admissible* invariant (similarity) solutions 
of equations (1) to (3) will be considered herein. Therefore, it is necessary to seek 
for and determine the form of such solutions. This can be done (see Theorem 3 of 
Ref. 3) by considering the invariance of equations (1) to (3) under the continuous 
one-parameter group given by f 


with m, q, r, s=0 and numerical parameter a. Now a complete set of absolute in- 
variants of the group (4) is given by 
6, and RM . . (5) 


Also, upon substitution of equations (4) into the system given by equations (1) to (3), 
it is found that these latter are constant conformally invariant under the group (4) if 


q=-m,r=—-m and J=-2m.. 7 (6) 


Hence, by the Theorem 3 of Ref. 3 and the use of equations (5) and (6), the sub- 
stitutionst 


py? 


PO . ; (7) 


v v 
4 = P= 
u= (4), RY: P 


will reduce equations (1) to (3) to a system of ordinary differenuial equations. The 


*Determined by the applicable boundary conditions. 

{This is a generalisation of the group of uniform expansions. 

tThese are the functional forms assumed by Yatseyev'!), equivalent to Squire’s) assumed 
functional form for the stream function. 


August 1956 229 


in 
(21) 
OrA 
in 
x Z 
of 4 
air 
m- | 
ted 
ing 
(1) 
(3) 
ity 
re 
rly 


A. J. A. MORGAN 


functions U (6), V (6), and P(@) then define the invariant (similarity) solutions of 
equations (1) to (3). In equations (7) the factors v and pv* have been introduced so 
as to make the functions of 6 dimensionless. In fact, U (@) and V (@) may be inter- 
preted as defining the local Reynolds numbers of the radial and tangential velocity 
components, and may then be referred to as Reynolds velocities. 


The resulting system of ordinary differential equations is 
U’V -U*?—V*= -2P+U”+U’ cot 6—2U —2V’ - 2V cot 6, . (8) 
VW=P'+V"+V’ cot 0+2U’ V cosec? 4, . (9) 


where the primes denote differentiation with respect to 9. 


Before proceeding further, that class of boundary conditions which can be cast 
solely in terms of the functions U (@) and V (4) must be isolated. If inner and outer 
conical boundaries are located at 0=A and @=A respectively, then the most general 
boundary conditions admissible under the restrictions (7) on the functional form of 
the velocity components are 


U(A)=U,, U(A)=U,, VO)=Vi, 
Hence, u and v will vary as R~' along the boundaries and be singular at R=0. 


Thus, suction or injection velocities, of the special form given by equations (7) 
and (11), at any two given angles with respect to inner and outer conical boundaries, 
may be prescribed. 


For example, the assignment of a positive number for V, is equivalent to pre- 
scribing a normal suction (i.e. fluid being removed from the space within the two 
conical boundaries) velocity at the outer conical boundary. 


There remains to be discussed the case where the inner boundary is removed and 
the fluid fills the whole space within the outer conical boundary. In this case two 
boundary conditions, U, and V,, can be assigned at, say, =A. On the axis, 6=0, 
it is evident that singularities in U (#) and V (@) cannot be admitted since this would 
imply the existence of a line source or sink within the fluid; hence, by symmetry. 
V (0)=0 and (dU /d6)s_,=0; U (0) can, however, be a finite positive or negative 
quantity. 


3. Solution of the Governing System of Ordinary Differential 
Equations 


On eliminating P (#) between equations (8) and (9), there results the expression 
(U’VY =(U”Y +(U’ cot +(U?Y +2U’, 


where the primes denote differentiation with respect to 6. By the use of the 
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continuity equation (10), equation (12) can be integrated three times in succession to 
yield a relation between U (4) and V (4), namely : —- 


U sin’ 6+2V sin 6 cos 6+4V? sin? 6= —4K sin?6+Lcosé6—M, . (13) 


| where K, L, and M are constants of integration. Again by the use of equation (10), 
(13) gives the following first order non-linear ordinary differential equation for 
V(@):— 


qo Sin’ @—V sin 6 cos 3V* sin’ 6=4K sin? 6—Lcos6+M. 


From equations (10) and (13) it can also be deduced that 


dv 


and 


The solution of equation (16), for arbitrary values of K, L, and M, can be 
| obtained by an application of the transformation 


V sin 9, (1 - and €=cos? 


t(g) dé 2 (17) 


The function 7 (€) is found to satisfy a hypergeometric differential equation so 
| that, if y; is not an integer™, 


Nij ({)=F (2i;, Bis; Vis €)+Cé! Yi F (3; | Vis (18) 


Where i, j=1, 2 and, because of the first of equations (17), only one constant of in- 
tegration, C, has been included in the expression. 


The four possible solutions ;;(€) given by equation (18) arise from the choice 
made for the roots of the two quadratic equations which it is found that p and q 
must satisfy. These roots are 


2p,=1+{1-4(L+M)}, 19) 


2q,=1+[1+4(L-M)}, 2q.=1-[1+4(L-M). (20) 


August 1956 231 


of 
so 
ter- 
city 
(8) 
(9) 
10) 
ast 
iter 
ral 
of 
|| 
11) 
(7) | 
es, 
wo 
nd 
WO 
ld 
ry. 
on 
he 
rly 


A. J. A. MORGAN 


If i, j=1, 2, then the parameters of the hypergeometric functions in equation (20) B qua 


are given as case 
con: 


The Reynolds velocities V (#) and U (@) are then, on using equations (10) and & s0 1 
(17), given by flov 


V (6)= 2p, tan — 2q; cot d (logen: sin (i,j=1.2), . (24) thr 
2 equ 
ma 
U (0)= 2p, ~ 24, +25 gett (23) 
5. 
dU (0) {logeni (©)} 


It should be noted that the right hand side of equation (24) always represents the F vis 
same function of @ irrespective of the choice made for p; and gq; since, as can be F 

shown from equations (14), (17), (19), and (20), the »,; (€) are inter-related by the Pr 
transformations Ui 


The same arguments apply to the right hand sides of equations (25) and (26). Hence, 


the choice of the pair (p;, qj) is purely a matter of convenience. The expressions 
(24) to (26) then complete the formal solution of the system of equations (8) to (10). 


4. The Mass Flow Rate 
The mass flow rate across a spherical surface intersected by the two bounding 

cones situated at 6=A and 6=A (A>A) is, with the aid of the continuity equation 
(10), given by 
wi 


A 
M =—2rpyR | (V sin 0) —2xpyR (V, sin A V; sin A). . (27) 


Thus, M, for the class of laminar viscous flows under consideration, is directly pro- 
portional to R and is in addition dependent on the magnitudes of the suction/ 
injection Reynolds velocities V, and V;. There will be a net mass outflow (i. 
away from the cone’s vertex) if M > 0 and a net inflow if M <0; that is, if the dl 
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quantity (V, sin A — V; sin A) is less or greater than zero respectively. In the special 


case where the boundaries are impermeable (V,=V;=0), M=0 and the flow could 


consist of a combination of out- and back-flow regions. 


For the case where the viscous fluid fills the space within one bounding cone 
situated at @= A, equation (27), since V (0)=0, is replaced by 


' so that there will be a net mass inflow if V, >0, a net outflow if V,<0, and no net 


flow if V,=0. A simple calculation will show that the value of M given by 
equations (27) and (28) is exactly equal to the mass rate at which fluid is introduced 
through the sides of the cones. The singularity in R of the velocity components (see 
equations (7)) at the origin is not sufficiently strong to contribute anything to this 
mass flow rate. 


5. A Class of Laminar Viscous Flows Within One Bounding Cone 


First, the following result can be shown to hold : — 


A necessary and sufficient condition for equation (14) to represent laminar 


| viscous flows within one bounding cone is that L=M= - K. 


Proof of the Sufficiency: Since on the cone’s axis (@=0), V (0)=0, dU/dé=0 and 
U (0) must be finite, an evaluation of equations (13) and (15) at @=0 shows that 
L=M and L=~-K respectively. This verifies the sufficiency of the condition 
L=M= -K. 


Proof of the Necessity: If L=M= - K, the application of the transformations 


_ 1 di 
V (6)= dé sin #, U (6)=2- =cos* 5, (29) 
will reduce equation (14) to the form 
4g? 


dé 


which, depending upon the value of K, has the following solutions'” 


C, cos (b log.é) + C, sin (b log.é), when b? = +1) > 0, or K < 1.(31a) 
, when b?=4(K+1)>0,orK >-1, (316) 
(C,+C, log. é ,whenK=-1. .  . 


For each of equations (31) it can be verified that at 6=0 (€=1), V (0)=0 and 
dU/d#=0. Hence, this proves the necessity of the condition L=M= - K. 
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ORDINATES 
y=V,ton A ,x =(k ton? -2U,] ton@a + 
y=V,tanaA ,x = -2U;,ton7A, A# 
y=Vitond # = 


Assignable boundary conditions. 


It should be noted that, in view of the first of equations (29), only one constant 
of integration C=C,C.,”' or C=C,~'C, (whichever is preferable in use) is available 
for satisfying one of the boundary conditions V (A)=V, or U (\)=U,,. 


On using the foregoing result, equation (15) reduces to 


dU 6 
tan; +V(U +2) 


From equation (13) it is found that K is uniquely determined by the relation 


A 
K tan? tan? A - 2U, tan \(V, tan A+4) . 63) 


A plot of this relation is shown in Fig. 2 and it is seen that there exist two possible 
values of V,, tan A for each given value of the quantity 


A 
(K tan? 2U,) tan? A, AA 5 
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0-6 
U(8,A) 
U(O,A) 


0-4 


Fig. 3(a). 


Radial Reynolds-velocity component of viscous flow within cones of semi-vertex angle .\, 
injection V (.\, .\) as in Fig. 3 (b). 


As an example consider the situation when K=0 and U,=0. From equations 
(29) and (315) it is seen that 


sin 6 
6) — 
V (0) and U (4) 


Therefore, U (A) will vanish if 


\ 
‘ 


C=00, or C= 


The first case (C= 00) yields that branch of the solution which gives the identically 
quiescent flow V (6)=U (6)=0. Thus, this is the only flow of the form (7) which 
will satisfy no-slip conditions at the conical boundary. The second value of C gives 
that branch of the solution where V,,= - 4 cot A (as can be verified from Fig. 2) and 
U (0)=2 cos A cosec* (.\/4). The values of U (6) and V (@) for this latter flow are 
shown in Fig. 3. 


In general, therefore, the two values of C required to satisfy the condition 
U (A)=U, correspond, one each, to the two possible flows associated with each value 
of K, A, and U,. 
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Fig. 3(d). 


Tangential Reynolds-velocity component of viscous flow within cones of semi-vertex angle .\, 
injection V (.\, A) as in the table above. 


0-4 


A 


6. A Class of Laminar Viscous Flows Bounded by Two Co-Axial 
Cones with a Common Apex 


For arbitrarily assigned values of U,, V,, U;, and V;, equation (13), when 
evaluated at =A and at =<, yields two simultaneous equations from which the 
constants of integration L and M can be obtained in terms of K and the parameters 
of the problem. These results can be substituted into equation (15) to give the 
slopes of the U (@)-profiles at the conical boundaries. 
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The relations for L and M in terms of K are used to express the pairs, say, 
(p,, 4,) and (z,,, 8,,) wholly in terms of K and the parameters U,, V,, Ui, Vi, A and 
A of the problem. The Reynolds velocities U (@) and V (6) are then obtained by the 
use of equations (24) and (25). The constants K and C must still be evaluated by 
the use of the two conditions V (A)=V;, V(A)=V,. Since, in the general case, the 
constant K occurs in the parameters (2, 8, y) of the hypergeometric functions enter- 
ing into the problem, it is seen that the determination of K and C will involve rather 
extensive computations. It may then be simpler to solve the original ordinary 
differential equation (14) by use of a differential analyser or an analogue computer. 


The situation just described is somewhat simplified when consideration is given 
to that class of flows within two bounding cones where the Reynolds velocities at the 
boundaries are prescribed in a manner such that 


(2U, tan?A + V, tan.\ (V, tanA + 4)] cos? A=0,. (36) 

and (2U; + V; tan A (V; tan A + 4)] cos?A=0, ‘ G67) 

where A >A and which includes the case U,=V,=U;=V;=0. If A or AH*7/2, 
then the expressions within the square brackets in equations (36) and (37) are equal 
to zero, and these relations are shown in Fig. 2 as Curves 2 and 3 respectively. It is 


seen that for each value of 2U;tan?A or 2U, tan? A<4 there exist two possible 
values of V (@) at the boundaries. If A=7/2 or A\=7/2, then 


These relations are depicted as Curves 4 and 5 in Fig. 2. 


When the conditions (36) and (37) are satisfied, 
L=-—4K (cosA +cos \), M= 4K (1 +cosA cos A), 


and, with the transformation (29), equation (16) becomes 


d*) 
S 


A A 
where 5 » C08" 


The two linearly independent solutions of this equation are also given in terms of 
hypergeometric functions and numerical methods must again be used for the deter- 
mination of the two remaining constants of integration K and C. 


Of particular interest is the case when the two conical boundaries are imper- 
meable and the no-slip conditions are imposed thereon; that is, U,=V,=U;=Vi=0. 
A simple physical argument will show that the only solution of equation (40) which 
will satisfy these homogeneous boundary conditions is one which gives an identically 
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quiescent flow: V (#)=U (6)=0. From equation (27) it is seen thgt there can be no 
net mass flow across that portion of any spherical surface, centred at the origin, with- 
in (and intersected by) the two bounding cones. As previously mentioned, the 
singularity in R at the origin of u and v (see equations (7)) is not sufficiently strong 
to produce any mass flow into the space bounded by the two cones. Therefore, there 
can be no forces acting on the fluid within the two bounding cones; hence, the fluid 
can only remain perfectly still and the flow must be identically quiescent. Such being 
the case, the only possible value of K which will satisfy these conditions, for 
arbitrary A and A, is K=0. 


7. Discussion 


It might be thought that the class of viscous flows of equation (7) is the axially- 
symmetric physical analogue of the two-dimensional viscous (“vergent”) flows 
between non-parallel plane walls considered by Jeffery” and independently by 
Hamel". This analogy is suggested by the fact that in both cases the component 
velocities are inversely proportional to the distance from the origin. However, if 
this analogy were to hold, the solutions for the axially-symmetric case would exhibit 
properties similar to those encountered for vergent flows if, in both cases, the no- 
slip conditions are satisfied at the boundaries. It has been shown that the only 
axially-symmetric viscous flows of the class of equation (7) satisfying these boundary 
conditions for one or two bounding cones are those which are identically quiescent. 
Thus, flows of the form of equation (7) are not physically analogous to “ vergent” 
planar flows. This can also be seen in a more fundamental fashion. In the vergent 
flow case the mass flow rate across any arc of a cylindrical surface centred at the 
origin is independent of its radius. For the class of viscous flows of equation (7) the 
mass flow rate across any spherical sector centred at the origin is proportional to its 
radius. This clearly indicates that these two flows cannot possibly have similar 
properties. 


In view of the preceding work the problem of a “ jet emerging from a hole in a 
plane wall” considered by Squire™ can now be given a different interpretation. The 
flows considered by Squire® are those which would be produced by injecting fluid 
parallel to the wall and towards the origin such that, in the present notation, 
U=U, and V=V,=0 at 6=A=7/2. The analysis of Section 5 is applicable and, 
therefore, by equation (33), K =2U,, where U, is negative. 


If slip-type flows are excluded, it is difficult to visualise the means by which 
fluid could be injected parallel to conical, or for that matter any other, boundaries. 
However, the solutions presented in the body of the paper may be used to determine 
the flows within one or two bounding cones in those cases occurring in practice 
where injection velocities of a form approximated by equations (7) and (11) are 
encountered or specified. 


The present work leads to an interesting negative conclusion. The result of the 
investigations for the case of one or two impermeable bounding cones indicates that 
none of the solutions of the Navier-Stokes equations (1) and (2) invariant under the 
group (4) can be used to describe incompressible laminar viscous flows within 
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conical diffusers; ‘“‘ None” is asserted because equations (7) are the only possible 
form for these solutions. Further investigation of Stezkin’s“® formulation might 
give useful results in this direction. However, the present work does indicate that 
a better understanding of the flow within conical diffusers could be experimentally 
obtained by measuring both the radial (wu) and tangential (7) velocity components of 
such flows. To the author’s knowledge, only measurements of the radial velocity 
component have been reported upon in the literature. 
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Note on the Wave Drags of Untapered Wings 
with Streamwise Tips at Zero Lift 


A. D. YOUNG 


(Queen Mary College) 


SumMMaARY: Detailed calculations based on linearised supersonic flow theory 

have demonstrated in a number of cases the equivalence between the wave 

drags of yawed or swept wings of constant section and zero lift and the 

corresponding values in two-dimensional flow. These results and others are 

shown to follow readily from a general theorem for which a simple proof is 
given. Some interesting deductions are illustrated in Fig. 3. 


1. Introduction 


It has long been established that the wave drag at zero lift of a rectangular 
wing of constant symmetrical section is, according to linearised theory, the same 
as in two-dimensional flow provided that the Mach cone from one forward corner 
does not cut the other streamwise edge, i.e. provided AB>1, where A is the aspect 
ratio and B=/(M,?—1). Again, for a semi-infinite sheared wing of constant 
symmetrical section with rear streamwise tip it is known that the net effect of the 
tip on the drag is zero whether the leading edge is subsonic or supersonic, i.e. the 
drag is the same as if the wing were part of an infinite yawed wing, while the same 
is true for a semi-infinite sheared wing with forward streamwise tip provided that the 
leading edge is supersonic. For sheared wings of finite span and constant section 
and for swept wings (i.e. of plan form) of constant section, similar results in 
certain cases have been found as a result of detailed and sometimes involved 
calculations*. These results and others can, however, all be deduced from one 
general theorem which is very simply proved and which, as far as the author is 
aware, has not been previously demonstrated. 


*It was one such calculation by G. Haynes at the College of Aeronautics that stimulated 
this investigation. 


Received July 1955. 
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WAVE DRAG OF SWEPT WINGS 


NOTATION 
aspect ratio 
V(M,?—1) 
free stream Mach number 
angle of sweep of leading edge 
Mach angle 
rectangular Cartesian co-ordinates 
see Fig. 2(b) 
free stream velocity 
ordinate of wing section upper surface 
wing section chord 

#,0+80 angles defining triangular strip of wing (see Fig. 2(b)) 


S(x,) = Jz dx 


Bcot A. 


2. Statement of the Theorem and its Proof 


Figure 1 illustrates the plan forms for three cases of semi-infinite wing of 
assumed constant symmetrical section with streamwise tip at zero incidence. In 
(a) the leading edge is supersonic and the tip is a forward one, in (b) the leading 
edge is supersonic and the tip is a rear one, and in (c) the leading edge is subsonic 
and the tip is a rear one. In each case it is visualised that a semi-infinite extension 
is brought up to the wing with the same sections in the direction of flight as the wing 
and with a streamwise tip that is made coincident with that of the wing, but the 
leading edge angle of wing and extension may in general differ. 


LIMITING 
DIRECTIONS 
FOR EXTENSION : 
LEADING EDGE LIMITING 
TIONS 
FOR EXTENSION 
LEADING EDGE LEADING EDGE 


(¢) 


Fig. 1. 
Sketch illustrating cases for which the theorem applies, 
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EXTENSION 


Fig. 2. 
Illustration of wing and extension. 


Theorem: In the cases illustrated in Fig. 1 the net effect on the drag of the 
wing due to the presence of the extension is zero, provided that the leading edge of 
the extension is supersonic if it has a rear vertex, while its leading edge can be 
either supersonic or subsonic if it has a forward vertex (the limiting directions are 
illustrated in Fig. 1). 


Proof: The theorem will be proved for the case of the wing with rear vertex. 
but the argument applies unchanged to all the other cases. 


Figure 2(a) illustrates the wing and extension. The wing section is assumed 
to have pointed leading and trailing edges. The effect of the extension on the flow 
over the wing is clearly confined to the region OAB where the angle AOB is », the 
Mach angle. The extension can be regarded'':*) as equivalent to a system of semi- 
infinite source planes, each of constant strength per unit area and bounded by the 
x-axis along OA and by a line parallel to the extension leading edge. Thus, the 
particular source plane bounded by the line parallel to the extension leading edge 
and passing through x=.x,, y=z=0, and associated with the narrow strip defined by 
similar parallel lines passing through the points x=x,+Ax/2, y=z=0, say, is of 
strength 2V, (0°z,/0x?),..,Ax where V, is the undisturbed stream velocity, and z 
is the ordinate of the upper surface of the wing section, while in addition there is a 
source plane beginning at the extension leading edge of strength 2V, (0z,/0x),-, and 
another beginning at the extension trailing edge of strength —2V, (0z,/0x),-., where 
c is the section chord. It can be checked that this system of source planes is 
equivalent to the extension, since the net source strength per unit span and unit 
distance parallel to the x-axis at any point is then 2V, (0z,/@x); as required by 
linearised theory. 
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WAVE DRAG OF SWEPT WINGS 


Each of these source planes will produce a conical flow field over the wing with 
vertex at the junction of the source plane leading edge and x-axis, and the total 
effect of the extension on the flow over the wing will be the sum of these conical 
flow fields. 


Consider now the source plane bounded by the line through x=x,, y=z=0, 
where 0< x, <c, and consider its influence over a narrow section of the wing, 
triangular in plan form and bounded by the lines through (x,, 0, 0) making angles 
6 and 6+60 with the x-axis, where 50 is small (see Fig. 2(b)). The contribution of 
the conical field due to the source plane to the pressure distribution over the wing 
will be a function of 8, independent of distance from the point (x,, 0, 0) and, for 
small enough 46, it may be taken as constant over the triangular region considered. 
Hence its contribution to the drag of this triangular strip of wing will be of 
the form 


c 


OZ. 
constant x V, Ax g (0) (=) (x—x,) dx, 


where g (4) is a function of #, independent of x. Hence, since z(c)=0, this contri- 
bution to the drag may be written 


constant x V, Ax g (9) 56 [S (c) S(x,)] 


where | dx. 


0 


It follows that the total contribution of the extension to the drag of the wing 
can be written 


c 


constant x V, do | d+ (2%) _ 


0 0 0 


Where the last term derives from the source plane bounded by the extension leading 
edge of strength per unit area 2 (0z,/0x),-,. 


This expression for the total contribution immediately reduces to 


c 


2. 1 0z,? 02s 
constant x V, | — (=) S(c)+ dx + 
0 0 


since z, (0)= z, (c)=0. 


Hence the theorem follows. 
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3. Some Deductions 


Consider the cases illustrated in Fig. 1. If in each case the extension is made 
a straight continuation of the wing, the combination forms an infinite yawed wing. 
However, the theorem has shown that the drag of the semi-infinite wing is unaffected 
by the presence of the extension; it follows that its drag must therefore be the 
same as if it were part of an infinite wing. Where the leading edges are supersonic 
(Figs. l(a) and 1(b)) the flow is in any case that appropriate to an infinite yawed 
wing, except in the regions included in the Mach cones from the forward edge 
of the streamwise tip. It follows that for those regions the overall drag is the same 
as would be obtained if they were part of an infinite yawed wing, although the 
detailed pressure distribution is different. With the subsonic leading edge (Fig. 1(c) ) 
the wave drag of the infinite wing is zero and hence the total drag of the region 
included in the tip Mach cone is zero. 


Now consider a sheared wing with streamwise tips and supersonic leading edge 
of finite span sufficiently large for the Mach cone from one tip not to intersect the 
other tip, 


i.e. AB(1-1*)>1, 


where m=Bcot A, and A is the angle of sweep of the leading edge. Suppose a 
semi-infinite extension with the same leading edge sweep is brought up to the rear 
tip. Then the argument of the theorem applies unchanged and it is deduced that 
the extension makes no change to the drag of the wing. It follows that the overiull 
drag coefficient of the wing must be the same as that of the infinite sheared wing 
with the same leading edge sweep. The result for the rectangular wing is a particular 
case of this more general result. 


If the sheared wing of finite span has a subsonic leading edge then, provided 
that the Mach cone from the rear tip does not intersect the forward tip [i.c. 
AB (1+1/m)=> 1], the same argument can be applied to produce the conclusion 
that the drag of the wing is the same as if it were part of a semi-infinite sheared 
wing with forward tip. 


Next consider the extension to be the image of the wing in the x-axis, so that 
the wing and extension form a swept or V wing. From the cases illustrated in 
Fig. 1 it will be clear that this can only be done within the validity of the theorem 
when the wing leading edge is supersonic. Suppose first that the wing is semi- 
infinite. Then it follows from the theorem that the drag of each half wing of the 
swept wing combination is unaffected by the presence of the other, and hence the 
total drag is just twice the drag of the semi-infinite sheared wing. Since the latter 
has been shown to be the same as if it were part of an infinite sheared wing, it 
follows that the drag coefficient of the swept wing combination is the same as that of 
the infinite sheared wing. 
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WAVE DRAG OF SWEPT WINGS 


(4) (b) (c) 
SEMI-INFINITE SHEARED WINGS 


(a) and (b) Supersonic leading edge. Drag coefficient same as for infinite sheared wing. 
(c) Subsonic leading edge. Drag coefficient is zero. 


\ 
\ 


@) 


FINITE SHEARED WINGS 
(d) Supersonic leading edge. Drag coefficient same as for infinite sheared wing. 


(e) Subsonic leading edge. Drag same as if wing were part of semi-infinite sheared wing with 
forward tip. 


AB 22 aB(I-m) 


(ink 


(f) (9) 
SWEPT WINGS 
(f) and (g) Supersonic leading edge. Drag coefficients same as for infinite sheared wing. 
Fig. 3. 
Summary of deductions from the general theorem. 
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The argument is readily extended to produce the same result for swept wing 
combinations of finite aspect ratio provided that the Mach cones from the leading 
edge of the centre line do not intersect either of the wing tips and vice versa, i.e. 
provided that A (1 - 1/m) >2/B. 


It is clear that the argument could be extended, if required, to wings made up 
of a number of spanwise sections of different sweep angles. 


Figure 3 summarises these deductions. 
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Review 


The Theory of Hydrodynamic Stability, C. C. Lin. Cambridge University Press, 1955. 
155 pp. Diagrams. 22s. 6d. 


This volume by Professor Lin forms one of the series of Cambridge Monographs 
on Mechanics and Applied Mathematics edited by Dr. Batchelor and Professor Bondi. 
It gives an admirable summary of the present state of our knowledge of this most 
difficult and complex topic in fluid dynamics. The author has himself made notable 
contributions to this field and is admirably qualified to expound the various theories 
which have been propounded. After an introductory chapter explaining the general 
method of approach, the author deals with plane Couette motion with special reference 
to the classical experiments of Sir Geoffrey Taylor on the flow between rotating cylinders. 
The next topic for consideration is the stability of plan Poiseuille motion in which 
the basic flow is a parabolic velocity distribution in a circular channel. The delicate 
matter of the relations between the different branches of the asymptotic as well as of the 
fundamental equation is tackled in this chapter by heuristic methods, a rigorous 
analytical approach being deferred until the concluding chapter of this book. 


Chapter 4 gives an account of the general theory of hydrodynamic stability and 
Chapter 5 is concerned with stability of the boundary layer on a flat plate. This is one 
of the most interesting topics for the aeronautical engineer in view of the well known 
experimental results obtained by Schlichting and Schubauer. A number of other topics 
are discussed in Chapters 6 and 7, including some stability problems of interest for 
astrophysics and geophysics. There is a copious bibliography and author index. There 
is no doubt that this will be the standard work on the subject of hydrodynamic stability 
for some time to come. The work covers not only the fundamental mathematical 
analysis but also a discussion of the physical mechanism of instability and the inter- 
pretation of experimental results in the light of theory. 


G. TEMPLE. 
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